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Preface 


The  purpose  of  this  thesis  was  to  develop  a  finite  element  model,  and 
from  there  develop  an  empirical  equation,  that  could  be  used  to  predict  a 
porous  materials  strength  based  on  the  average  pore  size  and  volumetric 
porosity.  Such  a  model  would  allow  a  material  designer  to  evaluate  how 
strong  a  porous  material  would  be,  before  producing  a  sample  and  doing 
extensive  testing. 

I  would  like  to  thank  Captain  H.  Gans  who’s  assistance,  especially  in 
learning  the  ins  and  outs  of  the  computer  code,  was  indispensable.  I  am 
indebted  to  Dr.  A.  Palazotto  for  his  suggestions  and  encouragement  to  pursue 
this  topic.  I  would  also  like  to  thank  Dr  R.  Ruh  of  the  Air  Force  Materials 
Laboratory  for  suggesting  pertinent  references  when  I  was  stuck  and  needed 
help. 


Paul  R.  Woodmansee 
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Abstract 

j 

Porosity  in  materials  greatly  reduces  the  strength  and  load  carrying 
characteristics  of  the  material.  Porosity  is  unavoidable  in  some  materials 
(particularly  ceramics)  and  is  sometimes  desirable  for  other  reasons,  such  as 
radar  reflection  properties,  diffusing  a  fluid  through  a  material,  and  adjusting 
the  heat  transfer  properties.  In  the  past,  the  effects  of  porosity  on  material 
properties,  particularly  a  materials  modulus  of  elasticity  (Young’s  modulus), 
have  been  determined,  by  fitting  test  data  to  one  of  three  theoretical  equation 
forms  (linear,  empirical  exponential,  and  Hassin’s  semi-empirical  equation, 
whichever  fit  the  data  best).  These  analytical  equations  all  require  that  the 
effects  of  porosity  on  material  properties  be  determined  experimentally  for 
several  cases,  and  then  other  cases  can  be  extrapolated.  This  limits  accurate 
prediction  before  a  material  can  be  fully  analyzed  and  produced. 

In  this  thesis  a  finite  element  model  using  MSC/NASTRAN  is  developed 
that  can  numerically  determine  the  material’s  modulus  of  elasticity  using  the 
limited  information  from  one  material  sample.  The  model  is  three  dimensional, 
and  simulates  pores  by  placing  small  elements  that  are  non-load  bearing  into 
the  structure.  These  voids  are  randomly,  and  unevenly,  distributed  (using  a 
Poisson  distribution)  to  better  simulate  the  response  of  a  real  porous  material 
to  a  load.  How  much  this  porous  model  deforms  can  be  used  directly  to 
calculate  the  porous  modulus  of  elasticity.  The  model  is  shown  to  be 
repeatedly  accurate. 

The  same  finite  element  model  is  used  to  demonstrate  the  effects 


on  material  behavior  of  changing  the  pore  size  and  distribution.  From  this 
information  a  trend  was  noticed  and  an  empirical  equation  developed  that 
predicted  a  materials  elastic  modulus  based  on  porosity  and  pore  size. 


xi 


I.  Introduction 


Porosity  exists  in  almost  all  materials  to  some  extent.  Therefore  the 
effects  of  porosity  should  be  of  concern  to  any  material  developed  for  a 
design.  Porosity  is  of  particular  interest  to  a  designer  who  plans  to  use 
ceramic  materials,  because  it  is  unavoidable  in  ceramics.  Sometimes 
porosity  is  desirable  for  other  than  structural  reasons  such  as  heat 
transfer  properties,  easier  production  (bringing  lower  cost),  radar 
reflection,  et  cetera. 

When  dealing  with  porous  materials,  there  are  two  structurally 
oriented  concerns:  (1)  the  onset  of  a  pore  growing  into  a  crack,  and  (2) 
the  macroscopic  effects  of  the  the  pores  on  material  strength.  This 
thesis  is  concerned  primarily  with  the  macroscopic  effects  of  porosity. 

The  growth  of  pores  into  cracks  will  be  discussed  only  briefly  for 
completeness.  The  modulus  of  elasticity,  also  called  Young’s  modulus,  is 
the  major  parameter  of  interest  when  looking  for  the  macroscopic  effects 
of  porosity.  Therefore,  this  thesis  will  deed  primarily  with  the  effects  of 
porosity  on  Young’s  modulus. 

Porosity  is  defined  as  the  fraction  of  volume  occupied  by  non  load 
carrying  voids  in  the  material.  Past  work  on  the  effect  of  porosity  on 
Young’s  modulus  (which  will  be  talked  about  in  detail  in  Chapter  Three) 
dealt  with  fitting  a  curve  to  actual  data  for  a  specific  material.  Thus  for 
every  material  one  had  to  make  dozens  of  samples,  each  with  different 
porosities,  and  then  try  and  fit  a  curve  to  the  data  points  to  show  how 
Young’s  modulus  varied  with  porosity.  While  this  method  has  proven  to 
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be  reasonably  accurate,  it  is  also  time  consuming,  requiring  extensive  lab 
work.  Also,  there  was  no  way  for  a  material  engineer  to  predict  the 
results  of  changing  an  untested  material's  porosity  without  actually 
producing  the  material  and  testing  it.  Furthermore,  previous  work  did 
not  adequately  predict  the  effects  of  changing  the  pore  sizes,  or 
distribution. 

This  thesis  is  an  attempt  to  create  a  finite  element  model  that 
simulates  a  porous  material.  With  this  finite  element  model  the  porosity, 
pore  size,  and  distribution  can  be  varied  to  see  the  effects  on  the  overall 
material.  The  finite  element  model  uses  a  common  "off  the  shelf" 
computer  code,  MSC/NASTRAN.  The  computer  model  used  in  this  thesis 
only  used  linear  elastic  deformations,  but  it  can  be  extended  into  the 
nonlinear  regime.  Even  though  plastic  deformation  and  microcracking  are 
significant  effects  on  the  microstructure,  it  will  be  apparent  that  a 
reasonably  accurate  model  of  the  macroscopic  effects  of  porosity  can  be 
made  using  only  linear  elastic  calculations. 


II.  Objective 


The  ultimate  objective  of  this  thesis  was  to  devise  a  method  to 
predict  the  effect  of  porosity  on  a  material’s  modulus  of  elasticity.  This 
method  was  applied  to  another  material  to  predict  the  effect  of  porosity. 
Finally,  the  analytical  prediction  was  checked  against  a  fully  tested 
material. 

Several  predictors  (i.e.  equations  showing  how  Young’s  Modulus 
varied  with  porosity)  existed  already,  but  each  relied  on  testing  the 
material  at  various  porosities,  and  fitting  them  to  a  curve.  None  of  the 
previously  used  equations,  (18:458-459  and  7:327-329)  (Chapter  Three) 
took  in  to  account  pore  size  or  distribution. 

The  first  step  in  coming  up  with  a  predictor  was  to  create  a  finite 
element  model  that  simulated  the  microstructure  of  a  porous  material. 
Previous  finite  element  models  of  porous  materials,  (11:909-913  and 
6:2-15),  used  two  dimensional  circular  voids  evenly  spaced  in  a  material. 
Instead,  a  three  dimensional  model,  with  cubic  voids,  in  a  random  and 
uneven  distribution  was  chosen  for  this  analysis.  Also,  unlike  previous 
works,  linear  elastic  deformation  was  chosen  for  this  analysis.  Cracking 
and  plastic  deformation  were  not  considered.  Fortunately  the  model 
selected,  after  some  refinement,  did  a  reasonably  good  job  of  predicting 
the  properties  of  the  studied  materials. 

The  model  was  then  used  to  look  at  various  aspects  of  porous 
material  behavior.  Stress  concentrations,  the  effects  of  changing  the 


pore  distribution,  the  effects  of  varying  the  pore  size,  and  an  equation 
to  estimate  these  effects  on  the  modulus  of  elasticity  were  all  examined. 

The  equations  derived  from  this  thesis  allow  an  engineer  to  get  an 
estimate  of  the  effects  that  porosity  will  have  on  a  material.  However 
while  the  results  are  reasonable,  there  is  still  the  possibility  of  error  in 
the  estimates.  Some  lab  testing  should  be  done  to  confirm  the  theoretical 


results. 


III.  Background 


This  section  discusses  some  of  the  previously  used  equations  for 
Young’s  modulus  as  a  function  of  porosity,  and  how  they  were  developed. 
Volumetric  porosity  is  determined  by  comparing  the  theoretical  material 
density  with  the  actual  total  density.  Also  discussed  in  this  chapter  is 
the  stresses  and  strains  acting  on  a  pore,  and  the  effects  of  void  growth 
and  coalescence. 

Elastic  Modulus  versus  Porosity  Equations 

Linear  Equation.  The  linear  relationship  between  porosity  and 
modulus  of  elasticity,  Young’s  modulus,  is  the  oldest  and  still  the  most 
commonly  used  expression.  The  following  equation  is  most  common  way  to 
express  this  linear  relation  (16:78-79  and  18:458-459); 

£  *  if  0(  1  -  aP)  (1) 

where 

E  =  Elastic  modulus  (Young’s  modulus) 

Eo  -  Elastic  modulus  without  any  porosity 
a  =  slope,  an  empirical  constant 
P  =  volume  fraction  porosity 

The  values  of  a  (slope)  and  Eo  (zero  porosity  modulus)  are 
determined  experimentally  by  taking  several  material  specimens,  each  of 
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different  porosity,  and  finding  the  elastic  modulus.  The  accuracy  of  this 
linear  equation  is  often  quite  good.  For  example:  in  Petrak  (16:78-79) 
eight  tests  of  porous  cobalt  oxide  (CoQ)  were  done.  After  doing  a  least 
squares  curve  fit  using  Equation  (1)  he  found  values  for  a  and  Eo  that 
fit  the  data  with  only  a  2.56%  error.  It  is  worth  noting  that  Petrak 
(16:78-79)  went  on  to  curve  fit  the  data  to  the  empirical  exponential,  and 
semi-empirical  equations  (discussed  below)  and  came  up  with  even  smaller 
errors,  but  returned  to  the  linear  equation  because  of  it’s  simplicity  and 
accuracy. 

Empirical  Exponential  Equation.  While  investigating  the  microstruc¬ 
ture  features  on  the  mechanical  properties  of  hot-pressed  ceramic  oxides, 
Spriggs  and  Vasilos  (21:187)  found  that  the  effects  of  porosity  on  the 
elastic  modulus  for  alumina  and  magnesia  fit  very  closely  to  an 
exponential  expression.  The  exponential  expression  was  of  a  form 
previously  suggested  by  Knudsen  (12:376-387)  and  by  Duckworth  (4:68). 
The  expression  is: 


(2) 

where 

E  =  Elastic  modulus  (Young’s  modulus) 

Eo  =  Elastic  modulus  without  any  porosity 
e  =  the  Napier  number,  2.71828... 
b  =  an  empirical  constant 
P  =  volume  fraction  porosity 
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Since  it  was  first  introduced,  the  empirical  exponential  expression 
has  been  used  extensively  by  people  trying  to  correlate  material  strength 
with  porosity.  For  some  materials  it  is  an  excellent  fit,  for  others  the 
error  is  more  significant  and  one  of  the  other  equations  must  be  used  to 
try  and  curve  fit  the  solution. 

Semi-empirical  Equation.  Unlike  the  linear  and  empirical  exponential 
equations,  which  were  primarily  guesses  of  a  curve  that  seem  to  fit  the 
data,  this  method  has  some  basis  in  a  derivation.  It  was  originally 
deduced  by  Hasselman  (10:452-453)  from  work  done  by  Hashin 
(9:143-150). 

Hashin,  (9:143-150),  wanted  a  theoretical  expression  for  the  change 
in  a  material’s  overall  elastic  modulus  if  it  were  really  a  two  material 
substance.  Hashin  assumed  that  on  the  microstructure  level,  a  primary 
material  made  up  most  of  the  volume,  and  a  secondary  material  occupied 
small  islands  in  this  sea  of  prime  material.  Next  Hashin  broke  the 
material  into  smaller  volumes  called  cells,  where  each  cell  had  exactly  one 
"island"  of  secondary  material  inside  of  it  (Figure  1).  Finally  Hashin 
assumed  that  each  cell  was  spherical,  and  the  secondary  material  inside 
of  it  was  also  spherical,  and  body  centered  (Figure  2).  This  last 
assumption  was  needed  to  allow  any  work  to  be  done.  However,  this 
assumption  is  untrue,  since  it  becomes  obvious  that  a  group  of  spheres 
will  not  fit  together  to  make  up  a  solid  material. 

Once  these  assumptions  are  made  the  problem  has  radial  symmetry. 
Hashin  then  derived  an  expression  for  the  modulus  of  elasticity  for 


7 


a  continuous  phase  containing  a  dispersed  phase  (the  "islands")  of 
spherical  particles,  where  each  phase  had  a  different,  but  known, 
modulus  of  elasticity.  The  equation  can  be  written  in  the  form; 


where 


1  + 


4 

>-h: 

P 

l-(A+  1) 

OH 

:»-y 

\p] 

(3) 


E  =  Elastic  modulus  (Young’s  modulus) 

Eo  -  Elastic  modulus  of  continuous  substance 

Ei  -  Elastic  modulus  of  dispersed  substance 

A  =  a  constant  (actual  value  dependent  on  many  factors) 

P  =  volume  fraction  of  the  dispersed  substance  to  total  volume 


Hasselman,  (10:453),  took  Equation  3  above  and  applied  it  to  porosity. 
If  it  is  assumed  that  the  dispersed  substance  is  a  pore  with  no  strength, 
compared  to  the  primary  material,  then  Ei  =  0  and  result  is  a  porous 
material.  The  equation  then  becomes: 


E 


i. 

i-(yi.  n?  j 


where 

E  =  Elastic  modulus  (Young’s  modulus) 

Eo  *  Elastic  modulus  without  any  porosity 
Asa  constant 
P  =  volume  porosity 


(4) 
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Since  then  common  usage  of  this  equation  has  been  changed,  because 
A  in  Equation  4  was  always  negative.  Letting  A  =  -  A  gives  the  more 
common  equation,  (16:78-79  and  7:327-329): 


£-  E , 


1  - 


AP 


where 

E  =  Elastic  modulus  (Young’s  modulus) 

Eo  =  Elastic  modulus  without  any  porosity 
A  =  a  constant 
P  =  volume  porosity 


(5) 


Equation  (5)  then  is  the  semi-empirical  equation  used  to  estimate  the 
effects  of  porosity  on  elastic  modulus.  The  constants  Eo  and  A  are 
determined  from  the  experimental  data.  This  can  be  done  most  easily  by 
using  the  following  matrix  equation  derived  by  Hasselman  (10:453)  to  best 
fit  data  to  this  curve: 


\ 


a  _  p  2 

j 


where 

n  =  number  of  data  points  taken 
Ei  =  Elastic  modulus  for  data  point  i 


z 


p, 

1-p, 


(6) 
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Pi  =  volume  porosity  for  data  point  i 

Eo  =  Elastic  modulus  without  any  porosity  (unknown  being  solved) 

A  =  a  constant  (unknown  being  solved) 

While  the  semi-empirical  Equation  5  is  difficult  to  deal  with,  it  does 
seem  to  be  more  accurate  than  the  linear  and  empirical  exponential 
relationships.  However,  once  again  no  provision  is  made  to  see  the 
effects  of  changing  pore  size  and/or  distribution. 

Empirical  Power  Equation.  For  completeness,  there  is  one  last 
equation  that  is  used  to  correlate  porosity  with  elastic  modulus.  In  a 
Soviet  paper  by  Bubnov  (2:89-92)  the  following  empirical  equation  derived 
by  M.  Yu.  Bal’shin  (1:166-172)  was  recommended: 

E  m  E 0(\  -  P)n  (7) 

where 

E  =  Elastic  modulus  (Young’s  modulus) 

Eo  =  Elastic  modulus  without  any  porosity 
m  =  a  constant 
P  =  volume  porosity 

Again  the  value  of  the  empirical  constants  m,  and  Eo  are  derived 
from  testing.  In  Bubnov  (2:89-92)  this  equation  was  found  to  closely 
correlate  to  the  tested  effects  of  porosity  on  the  elastic  characteristics  of 
Silicon  Nitride  (Si3  N4).  The  maximum  difference  between  the  predicted 
elastic  modulus  and  the  tested  value  of  elastic  modulus  was  less  than  3%. 
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Stresses  and  Strains  about  Pores 


Several  researchers  (mentioned  below)  have  done  finite  element 
analysis  of  pores.  They  were  primarily  concerned  with  modeling  one  pore 
and  the  stresses  and  strains  about  it,  or  sometimes  two  pores  and  there 
interaction. 

The  first  such  study  was  done  by  McClintock  (15:363-371)  in  which 
he  looked  at  an  isolated  void  in  an  infinite  region  of  non  hardening  rigid 
plastic  material.  He  assumed  an  infinite  cylindrical  pore  so  that  he  could 
do  a  two  dimensional  study.  McClintock  found  that  void  growth  rate 
for  this  case  depends  exponentially  on  the  mean  in  plane  normal  stress. 
Later  Rice  and  Tracy  (17:201)  treated  the  case  of  a  spherical  void  in  an 
infinite  region  of  non  hardening  rigid  plastic  materiaL  They  found  an 
exponential  dependence  of  void  growth  rate  on  the  mean  normal  stress. 

Gurson  (6:2-15)  wanted  to  account  for  void  interaction  in  an 
approximate  way,  so  the  voids  were  placed  in  a  finite  cell.  For  spherical 
voids  this  cell  was  a  cube,  and  for  infinite  cylindrical  voids  (2 
dimensional  analysis)  this  cell  was  a  square.  Gurson  then  modeled  the 
material  as  rigid  and  perfectly  plastic,  and  approximated  the  local  velocity 
field  as  a  linear  function  of  the  macroscopic  strain  rate.  Gurson  obtained 
some  extensive  results  for  the  yield  functions  of  a  porous  solid. 

More  recently,  using  modern  computational  methods,  Kitagawa  and 
Honke  (11:909-913)  did  a  finite  strain  numerical  analysis  to  investigate 
anisotropy  and  softening  due  to  void  growth.  They  used  a  two 
dimensional  circular  cell  model  much  like  Gurson’s  though  unlike  the 
Gurson  model,  it  was  not  isotropic. 


(a)  Model  A 


(b)  Model  B 


Figure  3 


Model 


Kitagawa's  Finite  Element  Models 


Figure  4.  Resulting  Pore  Deformations  and  Strains 


Figure  3  shows  the  finite  element  model  used,  and  Figure  4  shows  the 
strain  distributions  for  this  model  under  two  different  loading  conditions. 
The  results  from  model  A  show  that  the  area  surrounding  the  pore  does 
plastically  deform,  and  model  B  shows  large  stress  concentrations  and 
strain  deformation  between  two  pores. 

When  the  strain  deformation  between  two  pores  cause  the  stress  to 
reach  its  critical  yield  point,  the  material  will  elongate  and  the  two  pores 
will  coalesce  into  one  large  oval  shaped  pore.  If  the  pore  is  long 

enough,  then  it  can  be  considered  a  microcrack  and  crack  growth  theory 
applies.  Similarly,  if  an  isolated  pore  is  deformed  enough,  it  can  be  the 
beginning  of  a  microcrack. 
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IV.  Model  Creation 


Choosing  Elements 

Before  constructing  a  finite  element  model  of  a  porous  material,  it 
was  necessary  to  determine  what  solid  elements  would  give  the  most 
accurate  results  for  this  situation.  Therefore,  a  conversion  study  was 
carried  out  similar  to  a  small  version  of  the  final  model.  Forces  were 
applied  for  which  expected  displacements  could  be  calculated.  If  the 
actual  displacements  matched  the  theoretical,  then  the  elements  were 
assumed  to  be  correct. 

Figures  5  and  6  are  two  views  of  the  test  model  used.  The  model 
consists  of  six  cubic  elements  (MSC/NASTRAN  CHEXA  elements)  arranged  in 
a  cube,  called  a  cell,  with  a  smaller  cube  inside  of  it.  These  seven 
elements  make  up  a  superelement  which  is  copied  seven  times  to  make  up 
a  larger  cubic  model.  The  images  of  the  superelement  do  not  show  up  on 
Figure  5  and  6  which  show  the  entire  model. 

The  points  originally  on  the  Z=0  plane  are  constrained  to  the  Z=0 
plane,  but  can  have  X  and  Y  displacements,  i.e.  on  rollers.  An  upward  z 
force  is  applied  in  some  cases  to  the  corners  of  the  model,  and  in  other 
cases  as  a  pressure  to  the  entire  upper  model  surface.  The  whole  model 
then  displaces  (Figure  7). 

The  20  node  CHEXA  element  was  tried  in  the  model  first. 

MSC/NASTRAN  forced  the  use  of  3  point  Gaussian  integration. 
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Since  the  center  cube  of  the  cell  had  the  same  material  as  the  rest  of  the 
model  (i.e.  model  a  solid  cube)  one  can  easily  calculate  how  much  linear 
displacement  in  the  Z  direction  is  expected.  The  20  node  elements  did 
not  displace  enough;  they  were  too  stiff.  Next,  the  central  cube  was 
made  into  a  pore  by  giving  it  an  elastic  modulus  (E)  much  smaller  than 
the  material's.  This  caused  the  corners  of  a  cell  to  displace  too  much, 
and  the  sides  not  to  displace  at  all  (see  Figure  8).  In  short,  a  20  node 
CHEXA  element  was  not  adequate. 

Next,  eight  node  CHEXA  elements  were  tried.  For  the  solid  cube  case 
the  computer  results  matched  the  theoretical  results  exactly  for  both  the 
corner  forces  and  surface  forces  case.  Also,  when  the  centred  cube  of 
each  cell  was  made  into  a  pore,  by  eliminating  it  and  creating  a  void,  the 
model  deformed  a  bit  more  than  the  solid  block  case,  as  expected.  The 
eight  node  model  gave  the  identical  results  whether  one,  two,  or  three 
point  Gaussian  integration  was  used  (specified  with  the  PSOLID  card). 
Therefore  the  default  one  point  integration  was  used  in  the  final  model. 


The  Model 

The  model  used  in  the  main  part  of  the  analysis  is  similar  in  many 
ways  to  the  test  model  used  above.  The  overall  model  is  a  block,  and 
the  size  varied  with  each  test  case  (Appendix  B).  This  block  was  divided 
into  80  cubes  called  cells.  The  block  was  four  cells  by  four  cells  by  five 
cells.  More  cells  were  put  in  the  Z  direction,  as  that  was  the  direction 
force  would  be  applied.  Like  the  test  case,  the  model  was  fixed  to  the 
Z=0  plane,  and  a  distributed  force  applied  to  the  top  plane  of  the  model. 


19 


21 


Each  cell  was  then  made  by  an  image  superelement  of  one  of  the  primary 
superelements.  The  primary  superelements  were  a  collection  of  elements 
simulating  zero,  one,  two,  or  more  pores  inside  of  a  cell.  Figure  9  shows 
the  overall  model  with  only  the  primary  superelements  shown.  The  image 
superelements  are  placed  in  the  cubes  created  by  each  of  the  residual 
data  points  shown  in  Figure  9. 

Unlike  previous  models  which  assumed  all  pores  to  be  round,  this 
model  assumed  all  pores  to  be  cubes.  This  is  done  primarily  to  simplify 
the  model.  However,  Figure  10  shows  a  picture  of  a  porous  material 
(cobalt  oxide)  and  it  is  evident  that  the  pores  (black  areas)  are  irregular 
angular  shapes  situated  between  grain  boundaries.  Therefore  an 
argument  can  be  made  that  pores  are  more  accurately  (or  at  least  as 
accurately)  simulated  by  cubes  rather  than  by  spheres. 

Superelements.  To  simplify  processing  time,  this  model  makes 
extensive  use  of  superelements.  A  superelement  is  a  collection  of 
elements  that  make  up  some  substructure  of  the  whole  model.  In  this 
case  the  substructure  is  the  cell,  so  the  whole  model  is  partitioned  into 
a  number  of  cells. 

Superelements  are  processed  in  the  following  way  by  MSC/NASTRAN. 
Use  Figure  11,  based  on  material  from  the  MSC/NASTRAN  Handbook  for 
Superelement  Analysis  (5:1, 1-2, 1.1-3),  to  follow  the  steps  taken  to  process 
a  superelement  model. 

1)  First  the  overall  model  is  partitioned  into  the 
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Figure  9. 


The  Final  Computer  Model  Used  (Ref  20) 
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Figure  10.  Cobalt  Oxide,  x400  Magnification 


Given  the  loads  P^ ,  find  the  resulting  displacements  u^.  Assume  all  springs  have  unit  stiffness- 


C.  Stiffness  Seduction 


0.  Load  Reduction 

Pi-0.5  P2-l  P^-0.5 


©-* w-^ 

K2  similarly 


Pj-1.5  P4-3  P5-I.5 


Figure  11.  Steps  Followed  in  Superelement  Analysis  (Ref  5) 


24 


Assembly  and  Solution 


K  «  ♦  K2  ■  1-0 


u2“*  U3“*,  U4»2 

enforced 


u^«0  U2*2.5  113*4  u4*3.5  uj»0 


Total  Solution 


Figure  11.  Steps  Followed  in  Superelement  Analysis  (cont.  Ref  5) 
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various  superelements  and  residual  data  points,  step  B  in  Figure  11.  The 
residual  data  points  are  points  that  don’t  belong  to  any  superelement, 
but  are  used  to  connect  superelements  together. 

2)  All  the  points  and  elements  in  a  superelement  are  then  reduced  to 
one  equivalent  stiffness  matrix,  step  C  in  Figure  11.  In  this  reduction 
operation  each  superelement  is  treated  as  a  complete,  disjointed 
structure,  and  solved  independently  of  the  rest  of  the  structure  using 
standard  finite  element  matrix  reduction  techniques. 

3)  Any  load  on  the  superelements  is  then  reduced  to  either  a 
boundary  load,  or  an  equivalent  load  on  a  connecting  residual  data  point. 
Figure  11  step  D  shows  the  boundary  loads  Pi  and  Ps  being  solved  for. 

4)  Then  the  equivalent  stiffness  matrix  for  all  superelements  that 
connect  to  each  residual  data  point  are  summed  to  get  an  equivalent 
residual  structure,  step  E  in  Figure  11.  Any  load  applied  to  a  residual 
data  point  is  applied.  The  displacement  on  all  the  residual  data  points, 
and  thus  the  entire  residual  structure,  can  be  easily  solved  using 
conventional  finite  element  matrix  manipulations. 

5)  The  last  step  is  full  data  recovery,  step  F  in  Figure  11.  Here 
each  partial  solution  due  to  reducing  the  loads,  step  D,  is  added  to  the 
partial  solution  due  to  the  enforced  motion  of  the  residual  data  points, 
part  E.  This  produces  a  solution  for  the  motion  of  each  data  point  that 
is  identical  to  the  solution  produced  by  conventional  analysis,  show  in 
part  A  of  Figure  11. 

In  this  paper’s  model,  all  of  the  loads  and  boundary  restraints  are 
applied  to  residual  data  points,  so  steps  D  and  F  of  Figure  11  are 


superfluous.  Once  the  equivalent  stiffness  matrixes  for  each  superele¬ 
ment  are  found  (step  C)  they  can  immediately  be  summed  up  for  each 
residual  data  point  and  residual  data  point  displacements  found  (step  E). 

Image  superelements  are  superelements  that  are  repeated  several 
times  in  the  model.  This  reduces  calculations  because  once  the  primary 
superelement  has  been  solved,  each  of  the  image  superelements  can  use 
the  exact  same  equivalent  stiffness  matrix.  In  the  model  used  in  this 
thesis  there  are  80  cells  but  only  a  few  of  them  (9  or  less)  are  primary 
superelements,  the  rest  are  simply  copies  of  one  of  the  primaries.  This 
greatly  reduces  computer  run  time. 

Primary  Cell  Configurations.  Each  of  the  primary  cells  had  a  number 
of  pores  inside  of  it.  The  solid  cell  and  the  one  pore  cell  are  identical 
(Figure  12),  except  the  one  pore  cell  had  a  different  material  in  central 
cube  that  has  a  modulus  of  elasticity  one  million  times  less,  making  it  act 
like  a  non  load  bearing  pore. 

The  two  pore  cell  is  depicted  in  Figure  13.  The  two  pores  in  this 
cell  are  so  close  together  that  the  slab  between  them  will  possibly  be  in 
plastic  deformation  if  the  load  is  high.  However,  this  effect  was  ignored 
in  this  model  because  the  data  on  the  tested  material's  plastic  deformation 
was  unavailable.  The  final  results,  shown  later,  are  accurate  enough  that 
no  further  correction  was  needed. 

For  the  case  where  three  to  eight  pores  are  packed  into  a  cell,  it 
becomes  likely  that  the  pores  are  touching.  Figure  14  shows  a  cell  with 
eight  pore  sized  cubes  inside  of  it.  For  each  of  the  three  to  eight  pore 
cases,  a  number  of  the  eight  internal  cubes  equal  to  the  number 


Figure  13.  A  Double  Pore  Cell  (Ref  20) 


Figure  14.  A  Cell  with  Three  to  Eight  Pores  Inside  (Ref  20) 
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of  pores  were  assigned  a  material  with  almost  no  stiffness  to  simulate  a 
pore.  In  the  case  of  the  eight  pore  cell,  all  eight  central  cubes  have  no 
stiffness,  making  the  eight  pore  cell  a  cell  with  one  large  internal  pore. 

Porous  Cell  Distribution.  Now  that  cells  with  different  number  of 
pores  in  them  can  be  simulated,  how  should  the  model  be  populated?  For 
some  of  the  models  an  even  distribution  of  one  pore  per  cell  was  tried, 
but  for  most  cases  a  Poisson  distribution  (explained  below)  was  used. 

In  Haight’s  book  Applied  Probability,  "A  Poisson  distribution  is 
formed  by  normalizing  the  terms  in  the  exponential  sum  to  unity  (8:8).” 

To  do  this  one  takes  the  exponential  summation  formula: 


•*- 1  *x*\*2*  jr*3"”'  (8) 

where 

x  =  parameter  of  probability 
e  =  the  Napier  number,  2.71828... 

Then  divides  it  into  discrete  probabilities  which  must,  by  definition,  sum 
up  to  1.  The  discrete  probabilities  are  therefore  written: 


-x  1  2  -* 

,xe  ,  —  x  9  , 


3! 


3  -  x 

X  9  , 


(9) 


where  the  symbology  is  the  same  as  Equation  8  above. 


Usually  all  terms  in  a  probability  distribution  can  be  represented  by 


a  single  formula  with  a  dummy  variable  introduced  for  compact  notation. 
This  is  true  for  Poisson  distribution  which  can  be  written: 


xzo'x 

z\ 

where 

Z  ~  1,2,3,... 

x  =  parameter  of  probability 
e  =  the  Napier  number,  2.71828... 


(10) 


According  to  a  paper  by  R.D.  Thomson  and  J.W.  Hancock  which 
discussed  the  effects  of  non- homogeneous  void  distributions,  "The 
probability  that  any  cell  selected  at  random  will  contain  x  inclusions  is 
given  by  the  Poisson  distribution  (22:108-109)."  In  the  case  of  a  cell 
with  x  inclusions  in  it,  Thomson  and  Hancock  (22:108-109)  used  the 
following  Poisson  distribution  equation: 


P  (*) 


8*(X-  1)! 


where 

P(x)  =  probability  of  there  being  x  pores  in  the  cell 
x  =  number  of  inclusions  in  the  cell 
e  =  the  Napier  number,  2.71828... 

*  =  average  number  of  pores  per  cell 


(11) 
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The  average  number  of  inclusions  per  cell  is  obtained  from: 


-  v* 

*  V, 


where 

x  =  average  number  of  pores  per  cell 

Vt  =  volume  of  actual  cell  used 

Vx  =  volume  of  average  cell  with  one  pore 

And  the  volume  of  the  average  cell  with  one  pore  is: 


(12) 


(13) 


where 

Vx 

P 


=  volume  of  average  cell  with  one  pore 
=  pore  volume 
=  volume  porosity 


Therefore,  using  known  information;  pore  size,  cell  size,  and  volume 
porosity;  Equations  11,  12,  and  13  can  be  used  to  calculate  the 
probability  of  a  cell  having  x  number  of  pores.  Once  the  probability 
distribution  is  known,  it  is  a  simple  matter  to  allocate  a  certain  number 
of  cells  in  the  80  cell  model  to  match  the  probability  distribution. 
Appendix  B  gives  the  actual  distributions  used  for  each  case  tested. 
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Evaluating  Test  Results 


MSC/NASTRAN  solves  all  finite  element  models  in  terms  of 
displacements.  Therefore,  a  way  was  needed  to  turn  the  displacement 
back  into  a  value  for  the  modulus  of  elasticity.  Using  the  linear 
stress-strain  relationship: 


o  =  Ee  (14) 

or 

f-j  (IS) 

where 

a  =  linear  normal  stress 
c  =  linear  normal  strain 
E  =  modulus  of  elasticity 

The  definitions  of  normal  linear  stress  and  strain  are  given  by 
Equations  16  and  17  below: 


a  ■ 


f 

A 


(16) 
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(17) 
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where 


a  -  linear  normal  stress 
F  =  applied  force 

A  =  area  over  which  force  is  applied 
c  =  linear  normal  strain 
AZ  -  deformation 

L  =  original  length  before  deformation 

Putting  the  equations  for  stress  and  strain  (Equations  16  and  17) 
into  Equation  15,  and  applying  the  result  to  this  model,  you  can  get  the 
following  relationship. 


E 


FAZ 

AL 


Where 

E  =  calculated  modulus  of  elasticity 

F  =  total  applied  upward  force  to  the  top  of  the  model 
A  =  model  top  end  area 
AZ  =  displacement  of  the  top  of  the  model 
L  =  model  height 


(18) 


Equation  18  was  used  throughout  all  of  the  testing  to  calculate  what 
the  total  model's  modulus  of  elasticity  was  for  that  test.  For  the 
distributed  pore  cases,  since  the  cells  were  distributed  randomly,  the 
final  displacements  on  the  top  end  of  the  model  are  uneven.  The  actual 


top  end  displacements  are  given  in  Appendix  B.  To  calculate  of  modulus 
of  elasticity,  the  average  displacements  were  used  in  the  distributed  pore 


cases. 

Global-Local  Analysis 

In  a  complex  structure  there  are  areas  where  a  fine  mesh  of  finite 
elements  are  needed  to  give  detailed  stress  and  strain  information,  such 
as  bolt  holes,  corners,  etc.  Such  a  fine  mesh  is  unnecessary  for  the 
entire  structure,  which  can  use  a  course  mesh.  Therefore,  fine  mesh 
area's  are  created  in  local  areas  of  interest,  and  the  results  of  the 
course  models  are  applied  to  the  fine  mesh  as  boundary  conditions  to 
obtain  local  effects.  This  is  called  global- local  analysis  (14:1). 

Normally,  the  global-local  method  treats  the  entire  structure  using 
the  Rayleigh-Ritz  method  but  adds  finite  elements  in  regions  of  interest. 
The  Rayleigh-Ritz  functions  need  not  satisfy  as  many  boundary  conditions 
then,  because  there  are  more  elements  to  smooth  things  out  in  the  region 
of  interest.  Global-local  analysis  is  useful  in  this  context  when  the 
structure  is  regular  enough  for  classical  Raleigh  Ritz  treatment,  and  only 
some  parts  need  more  detail,  thus  reducing  the  number  of  degrees  of 
freedom  for  a  given  level  of  accuracy. 

Global-local  analysis  is  used  in  this  thesis  in  a  slightly  different 
way.  The  computer  model  of  this  thesis  is  a  small  portion  of  a  large 
piece  of  material,  and  the  global  local  principle  is  used  to  deduce  that 
the  larger  material  will  behave  the  same  way  as  the  small  sample.  For 
example;  most  of  the  tests  use  30,  20  micrometer  cells,  making  the  entire 


model  have  dimensions  of  only  0.8  by  0.8  by  1.0  millimeters.  This  is  much 
smaller  than  even  the  smallest  real  test  specimen,  but  the  global  local 
principle  says  that  the  entire  structure,  which  is  made  up  the  same  way, 
will  deform  proportionately  the  same  as  the  small  computer  model.  This 
saves  a  great  amount  of  computer  time  because  now  the  entire  large 
structure  need  not  be  analyzed  in  the  fine  detailed  way  that  the  small 
piece  was  analyzed. 


V.  Results 


Stresses 

Before  looking  at  the  effects  of  porosity  on  the  elastic  modulus,  it  is 
worthwhile  to  examine  the  effects  that  pores  have  on  the  stress 
immediately  around  them.  To  do  this,  several  sample  cells  have  been 
taken  out  of  the  model  and  stress  contours  plotted  using  SCRD-IDEAS 
(20). 

Figure  15  shows  the  displacement  of  a  cell  with  a  single  pore  (like 
the  one  show  in  Figure  12).  Notice  that  the  displacement  is  uneven  due 
to  the  effects  of  neighboring  cells.  Figure  16  shows  the  exterior  stress 
contour  on  that  same  cell.  The  effects  of  a  small  central  pore  are 
minimal. 

Figure  17  shows  the  external  stress  contour  of  a  two  pore  cell, 
similar  to  the  one  in  Figure  14.  Notice  that  the  stress  contour  lines  are 
much  closer  together  indicating  a  higher  stress  gradient  across  the  cube. 
Figure  18  shows  the  stress  contour  internal  to  the  cell  along  the  surface 
of  the  two  pores.  The  stress  builds  up  and  is  largest  along  the  thin 
wedge  between  the  two  pores. 

Figure  19  shows  the  external  stress  contour  of  a  three  pore  cell.  In 
sharp  contrast  to  this  is  Figure  20  which  shows  the  interned  stress 
contours  of  that  same  cell.  The  three  pores  are  the  cubes  with  the  few 


stress  concentration  lines  on  them. 
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Figure  21  shows  the  external  stress  contour  of  a  four  pore  cell. 
Figure  22,  the  internal  stress  contours  of  the  same  four  pore  cell,  shows 
much  sharper  stress  gradients  than  Figure  21.  Thus,  the  stress 
gradients  close  to  the  pores  (Figure  22)  are  quite  high,  while  the  stress 
gradients  less  than  one  quarter  of  a  pore  radius  away  (Figure  21)  are 
much  less.  The  pores  in  Figure  22  are  the  four  cubes  with  few  stress 
lines  on  them,  while  the  solid  parts  have  a  tight  grid  of  stresses.  Figure 
23  and  24  show  the  external  and  internal  stress  contours,  respectively,  of 
a  five  pore  cell. 

The  trends  in  each  of  these  cases  is  clear  enough  with  these 
examples.  As  the  number  of  pores  rise,  both  the  stress  and  the  stress 
gradient  on  the  material  near  the  pore  rise  dramatically,  but  the  stress 
even  one  pore  diameter  away  changes  very  slowly,  if  at  all. 

Cobalt  Oxide 

If  the  porous  material  computer  model  was  going  to  be  useful,  the 
predicted  results  had  to  match  laboratory  tests  for  an  actual  material. 

The  material  initially  chosen  was  cobalt  oxide  (CoO)  because  of  the  wealth 
of  information  known  about  it  from  Petrak  (16:78-79).  Figure  10,  shown 
previously,  is  a  sample  of  cobalt  oxide  at  x200  magnification.  From  this 
picture  one  can  measure  the  average  pore  size  to  be  about  8  micrometers 
(i.e.  8  x  10®  meters).  This  average  ignores  places  where  two  pores  have 
combined  into  one,  since  the  model  takes  into  account  multiple  pores  that 
touch  each  other. 

First,  evenly  distributing  the  pores,  one  per  cell,  was  tried.  For 
this  case  the  toted  porosity  was  changed  by  changing  the  pore  size. 


Four  Pore  Cube  (Cell  49,  P0RE5-1) 


Max.  Stress  =  5.95 

Kez 

1  =  0.33  x  10-2  GPa 

2  =  1.27  x  10-2  GPa 

3  =  2.20  x  10-2  QPa 

4  =  3.14  x  10-2  GPa 

5  =  4.08  x  10-2  QPa 

6  =  5.01  x  10-2  QPa 


10-2  QPa 


Figure  21.  Four  Pore  Cell  Exterior  Stress  Contour  (Ref  20) 


Four  Pore  Cube  (Cell  49,  PORE5-1) 


Max.  Stress  *  5.28  x  10~2  GPa 


Key 

1  =  0.23  x  lO-2  GPa 

2  =  1.07  x  lO*2  GPa 

3  =  1.91  x  lO”2  GPa 

4  =  2.76  x  lO-2  QPa 

5  =  3.60  x  lO-2  GPa 

6  =  4.44  x  lO-2  GPa 


Figure  22.  Four  Pore  Cell  Interior  Stress  Contour  (Ref  20) 


Five  Pore  Cube  (Cell  45,  PORE5-1) 


Max.  Stress  =  7.49  x  10-2  GPa 


i 


1 

2 

3 

4 

5 

6 


Key 

=  0.96  x  lO-2  GPa 
=  2.05  x  10-2  GPa 
=  3.14  x  10-2  GPa 
=  4.22  x  10-2  QPa 
=  5.31  x  10-2  QPa 
=  6.40  x  10-2  QPa 


Figure  23. 


Five  Pore  Cell  Exterior  Stress  Contour  (Ref  20) 
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Five  Pore  Cube  (Cell  45,  P0RE5-1) 

Max.  Stress  =  5.19  x  10"2  GPa 

Key 

1  =  0.63  x  10“2  GPa 

2  =  1.39  x  10-2  GPa 

3  =  2.15  x  10-2  GPa 

4  =  2.91  x  10-2  GPa 

5  =  3.67  x  10-2  GPa 

6  =  4.43  x  10-2  GPa 


Figure  24.  Five  Pore  Cell  Interior  Stress  Contour  (Ref  20) 
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The  results  are  shown  on  Figure  25  (the  symmetric  distribution),  and  the 
raw  data  is  given  under  test  P0RE1-  in  Appendix  B.  Comparing  the 
results  of  this  case  to  the  actual  test  data,  also  shown  on  Figure  25,  one 
can  see  that  evenly  distributing  the  pores  causes  the  modulus  of 
elasticity  to  be  much  higher,  and  drop  much  slower  as  porosity  increases, 
than  the  real  test  values. 

Next  a  Poisson  distribution  (explained  previously)  was  tried.  This 
put  more  pores  in  some  cells  and  fewer  in  others,  while  keeping  the  pore 
size  of  8  micrometers  constant.  Figure  25  also  shows  how  this  turned 
out,  and  in  Appendix  B  the  data  is  given  under  the  P0RE2-  tests. 

Figure  25  shows  that  the  slope  using  distributed  pores  is  much  steeper 
than  the  case  for  evenly  distributed  pores  (i.e.  modulus  of  elasticity  for 
this  case  drops  much  faster  as  porosity  increases).  In  fact,  the  slope  is 
close  to  the  actual  test  results  once  they  are  linearly  curve  fitted,  but 
for  some  reason  all  the  computer  model  values  sure  high.  ThiB  is  unusual 
since  the  vsdue  of  zero  porosity  Young’s  modulus  (Eo)  used  for  the  solid 
elements  is  the  same  as  the  linear  approximations  intercept  point  for  zero 
porosity  (Eo  =  29.3  x  106  psia).  Thus,  one  would  expect  the  interception 
with  zero  porosity  to  be  about  the  same,  and  all  the  values  to  be  close  to 
the  same  since  the  slope  is  close. 

To  correct  this  problem,  the  value  for  the  modulus  of  elasticity  of 
the  solid  elements  was  dropped  to  26.1  x  10®  psia.  This  produced  the 
results  of  test  PORE2-  ;2  also  shown  by  the  shortest  dotted  line  on 
Figure  25.  The  results  are  close  to  the  actual  data. 
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0.00  5.00  10.00  15.00  20.00  25.00 

Porosity  (volumetric  percentage) 

Figure  25  Cololt  Oxide  Elastic  Modulus  vs  Porosity 


The  results  are  not  perfect,  in  fact  even  the  linear  least  squares  fit  line 
is  closer  to  the  real  test  data.  However,  these  results  were  computer 
simulated  using  only  the  value  of  zero  porosity  modulus  of  elasticity  (Eo) 
as  input. 

Silicon  Carbide 

Next  a  demonstration  that  the  model  worked  for  other  materials  was 
needed.  The  material  used  next  was  silicon  carbide  (SiC).  Ruh  (19:1369) 
states  that  silicon  carbide  has  a  reasonably  consistent  average  grain 
sizes  of  2.5  micrometers  (2.5  x  10**  meters),  and  gives  a  value  of  zero 
porosity  elastic  modulus  (Eo  =  440  Giga  Pascals)  as  the  linear  equations 
zero  porosity  intercept  point.  Figure  26,  and  PORE3-  cases  in  Appendix 
B,  shows  the  results  of  using  this  data  in  the  computer  model.  This  first 
run  was  slightly  high  just  like  the  first  cobalt  oxide  case,  but  lowering 
the  value  of  zero  porosity  elastic  modulus  (Eo)  to  425  Giga  Pascals  gives 
results  that  almost  exactly  matches  the  actual  test  data. 

This  demonstrates  that  the  model  was  accurate  for  more  that  one 
material.  It  also  shows  that  the  value  of  zero  porosity  elastic  modulus 
used  in  the  computer  model  should  be  about  5%  lower  than  the  linear 
equation  would  predict  it  to  be. 

Aluminum  Nitride 

It  has  always  been  difficult  to  get  a  close  correlation  between 
porosity  and  elastic  modulus  for  aluminum  nitride  because  aluminum 
nitride  pores  vary  greatly  in  size  and  shape  depending  on  the 
manufacturing  techniques.  Table  1,  from  Ruh  (19:1370),  shows  how  much 
pore  size  can  vary.  Notice  that  larger  pore  sizes  seem  to  cause 
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Porosity  (volumetric  percentage) 

Figure  26  Silicon  Carbide  Elastic  Modulus  vs  Porosity 


a  much  larger  elastic  modulus  even  if  the  porosity  is  about  the  same. 

For  example:  sample  6  has  a  larger  porosity  than  sample  5  and  a  larger 
(not  smaller  as  would  normally  be  expected)  elastic  modulus  due  to  much 
larger  pore  sizes. 


Sample 

no. 

Origin  * 

Porosity 

(%) 

Pore  size 
micro  m 

Elastic  Modulus 
million  psia 

1 

STK/-400/1700/5 

0.5 

2-10  (5) 

46.6 

2 

Fabricated  by  GE 

1.7 

1-3 

44.5 

3 

AT/ar/ 1 7 00/30 

4.1 

1-20  (10) 

40.8 

4 

AT/ar/1700/40 

4.7 

1-20  (10) 

42.3 

5 

STK/ar/1700/30 

5.4 

2-8 

38.6 

6 

AT/ar/2 1 00/60 

5.6 

5-50  (10) 

42.1 

7 

AT/ar/ 1750/ 15 

6.0 

1-5 

37.9 

8 

STK/ar/ 1700/20 

10.5 

2-6 

34.9 

9 

AT/ar/1850/30 

12.9 

2-10  (5) 

34.0 

10 

AT/+400/1600/15 

13.2 

5-120 

36.5 

11 

AT/+400/1700/15 

14.0 

5-120 

36.4 

12 

STK/ar/ 1600/ 30 

16.1 

1-10  (5) 

34.3 

13 

AT/ar/1950/30 

18.4 

2-20  (10) 

29.7 

14 

STK/ar/1600/30 

21.2 

2-20  (10) 

27.8 

15 

AT/+400/1700/20 

23.9 

3-250 

23.4 

*  =  Starting  powder/  Particle  size/  Hot-pressing  temp  (°C)/  Time  (min) 
STK  =  Stack,  AT  =  Atommergic,  ar  =  as  received,  GE  =  General  Electric 
Pore  sizes  in  parenthesis  used  as  actual  test  data  in  figure  27 


Table  1:  Microstructure  data  for  Aluminum  Nitride  (AIN) 


Aluminum  nitride  was  modeled  next  to  find  out  how  pore  size  affected 
the  macroscopic  material  behavior,  and  if  the  computer  model  simulated 
this  change  properly. 

First  cubic  pores  with  each  side  8  micrometers  long  were  tested, 
Appendix  B  test  results  P0RE4-  show  the  results,  using  a  zero  porosity 
elastic  modulus  (Eo)  of  46.7  million  psia.  The  results  fit  quite  close  to 
the  actual  test  data  given  in  Table  1  (Figure  27). 


V ) 

in 

(ft 

(ft 

0 

0 

0 

(ft 

(ft 

<D 

O 

U 

O 

O 

V 

L. 

<U 

i- 

0 

CL 

CL 

a 

0 

0 

a. 

Q. 

a. 

n 

O 

0) 

<u 

0 

V- 

O  0) 

*n  -*r* 

•O 

0) 

a; 

<U 

a) 

w  <u 

E 

E 

E 

<u 

<u 

test 
ro  m 

(/) 

E 

E 

0)  O 

0 

0 

JZ/  L_ 

O 

L- 

u 

v_ 

U 

O 

l- 

O 

—  0 
o*= 

0  E 

E 

E 

0 

to 

3  E 

3 

E 

E 

■4-J 

C  5 

0  0 

in 

0 

< 

00 

1 

m 

* 

* 

1 

1 

1 

1 

1 

l 

♦  /t  > 

'  s 

r 


/  j/ 

A/ 

,* 


//  J 
P  S' 


u/ 


0 

0 

O 

0 

O 

O 

0 

0 

0 

O 

0 

O 

q 

0 

d 

in 

d 

in 

d 

in 

d 

in 

rf- 

ro 

ro 

CnI 

CN 

Disd  uoi||i|Aj  ‘(3)  sn|npLM  oj^sog 


55 


0.00  5.00  10.00  15.00  20.00  25.00 

Porosity  (volumetric  percentage) 

Figure  27  Aluminum  Nitride  Elastic  Modulus  vs  Porosity 


Next  5  micrometer  pores  were  tested  (Appendix  B  test  results  PORE6- 
).  The  results  were  much  too  high  elastic  modulus  for  porosities  above 
two  or  three  percent.  The  reason  for  this  was  that  the  cell  size  was  too 
big  for  the  smaller  pore  size.  To  get  large  porosities,  i.e.  above  5 
percent,  large  conglomerations  of  up  two  ten  pores  in  one  cell  were 
needed,  and  few,  if  any,  cells  with  only  a  few  pores  in  them  were  used. 
Essentially  this  means  that  the  model  was  really  modeling  larger  pores, 
one  per  cell,  not  small  scattered  pores. 

This  was  solved  by  making  the  cell  smaller,  12  micrometers  instead  of 
20  micrometers.  This  made  distributions  of  zero,  one,  and  two  pore  cells 
much  more  prevalent,  and  eliminated  eight,  nine  and  10  pore  cells.  The 
results,  shown  in  Figure  27  and  in  Appendix  B  under  PORE7-  ,  were 
reasonably  close  to  the  actual  data  point  for  5  micrometer  pores,  shown 
on  Figure  27  by  an  asterisk.  Figure  27  also  shows  that  the  5  micrometer 
pore  results  are  almost  identical  to  the  8  micrometer  pore  results. 

If  the  10  micrometer  pore  results  were  not  significantly  different 
from  the  5  and  8  micrometer  pores,  then  the  model  wasn’t  properly 
accounting  for  pore  size.  The  10  micrometer  pore  tests  (P0RE5-  in 
Appendix  B)  shown  in  Figure  27  were  in  fact  higher  than  the  5  and  8 
micrometer  pore  results.  They  were  also  quite  close  to  the  actual  test 
results,  shown  by  stars  on  Figure  27.  The  10  micrometer  tests  showed 
that  the  model  did,  therefore,  show  the  effects  of  increasing  pore  size, 
i.e.  as  pore  size  grows  the  elastic  modulus  decreases  slower  as  porosity 
increases. 

Lastly,  15  micrometer  pores  were  tried  (PORE8-  in  Appendix  B)  to  see 


if  the  trend  continued.  Figure  27  shows  that  the  trend  not  only 
continues,  but  increases.  In  other  words;  as  the  pore  size  increases  the 
negative  slope  of  the  elastic  modulus  versus  porosity  line  decreases  at  a 
faster  and  faster  rate.  This  would  explain  why  there  was  very  little 
difference  between  the  5  and  8  micrometer  pore  results. 

Derived  Equation 

From  the  aluminum  nitride  model,  a  trend  can  be  seen  relating  the 
pore  size  with  the  rate  of  change  of  elastic  modulus  to  porosity. 
Therefore,  the  next  task  is  to  try  and  find  an  equation  that  describes 
this  trend.  This  was  done  by  taking  the  computer  generated  data  and 
trying  to  fit  it  to  the  previously  used  equations;  linear  (Equation  1), 
empirical  exponential  (Equation  2),  semi-empirical  (Equation  5),  and 
empirical  power  (Equation  7).  I  found  the  best  fit  to  be  with  Equation  5, 
the  semi-empirical  equation,  using  Equation  6  to  find  the  empirical 
constants.  Figures  28,  29,  30,  and  31  show  the  close  correlation  between 
the  computer  data  and  the  empirical  equation  for  pore  sizes  of  5,  8,  10, 
and  15  micrometers  respectively.  In  all  cases  except  the  10  micrometer 
pores,  the  error  is  less  than  1  percent.  Now  all  that  was  needed  was  a 
relationship  between  the  pore  sizes  and  the  semi-empirical  equation 
constants  to  get  a  good  predictor. 

To  find  this  relationship  the  semi-empirical  equation  was  rewritten  a 
bit.  Previously  the  semi-empirical  equation  was  given  as: 
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0.00  5.00  10.00  15.00  20.00  25.00 

Porosity  (volumetric  percentage) 

Figure  28  5  micro  meter  pore,  Al  N,  Elastic  Modulus  vs  Porosity 
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0.00  5.00  10.00  15.00  20.00  25.00 

Porosity  (volumetric  percentage) 

Figure  29  8  micro  meter  pore,  Ai  N,  Elastic  Modulus  vs  Porosity 
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Porosity  (volumetric  percentage) 

Figure  30  10  micro  meter  pore,  Al  N,  Elastic  Modulus  vs  Porosity 
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Figure  31  15  micro  meter  pore,  Al  N,  Elastic  Modulus  vs  Porosity 


where 


E  =  Elastic  modulus  (Young's  modulus) 

Eo  =  Elastic  modulus  without  any  porosity 
A  =  a  constant 
P  =  volume  porosity 


This  form  is  slightly  misleading  as  Eo  in  Equation  5  is  usually  not  the 
actual  zero  porosity  elastic  modulus.  So  another  parameter,  D,  was  added 
to  correct  this,  making  the  equation: 


E  •  DE , 


1  - 


AP 


l+(d-l)Pj 


where 

E  =  Elastic  modulus  (Young's  modulus) 

Eo  =  Elastic  modulus  without  any  porosity 
A  =  a  constant  parameter 
D  =  a  constant  parameter 
P  =  volume  porosity 


(19) 


In  Equation  19  the  value  of  Eo  is  the  actual  zero  porosity  elastic  modulus, 
and  the  parameter  D  modifies  it  according  to  the  pore  size. 


So  now  a  relationship  between  pore  size  and  the  parameters  A  and  D 
was  needed.  Table  2  shows  the  parameter  values  (calculated  using 
Equation  6)  and  pore  sizes  for  aluminum  nitride. 


Pore  size 

Calculated 

Parameter 

Parameter 

Micro  meters 

Eo 

D  * 

A 

5 

52.26 

1.1191 

3.7948 

8 

52.43 

1.1227 

3.7251 

10 

51.09 

1.0931 

3.12085 

15 

47.05 

1.0075 

1.51536 

*  =  Assuming 

actual  Eo  is 

46.7  million 

psia 

Table  2  Equation  19  parameters  for  aluminum  nitride 


Using  a  Shareware  program  called  CURVEFIT  by  Thomas  Cox  (3),  an 
equation  was  found  that  almost  perfectly  matched  the  relationship 
between  parameter  A  and  the  pore  size.  The  program  called  it  Hoerl’s 
Equation; 

>4  -  (P  1  )(P2i)(lM)  (20) 

where 

A  =  parameter  to  use  in  Equation  19 
L  =  pore  size  (length,  diameter,  etc.)  in  micrometers 
PI  =  a  constant,  in  this  case  0.8624 

P2  =  a  constant,  in  this  case  0.7461 

P3  -  a  constant,  in  this  case  1.831 

Figure  32  shows  how  nicely  the  Parameter  A  is  described  by 


Equation  20. 


#  Using  the  same  program  (3)  a  relationship  was  also  found  between 

the  parameter  D  and  the  pore  size.  A  Cauchy  distribution.  Equation  21, 
best  described  the  relationship  between  the  two  values. 


D  - 


1 


/*!(£  +  P2y  *  P3 


(21) 


where 

D 

L 

PI 

P2 

P3 


parameter  to  use  in  Equation  19 

pore  size  (length,  diameter,  etc.)  in  micrometers 

a  constant 

a  constant 

a  constant 


The  values  of  PI,  P2,  and  P3  for  Equation  21  that  best  fit  the  data 
in  table  2  are  shown  by  the  solid  line  in  Figure  33.  However,  If  the 
condition  that  D  goes  to  one  when  pore  size  goes  to  zero  is  imposed,  then 
the  constants  take  on  the  values  given  by  the  dotted  line  in  Figure  33. 
Neither  curve  exactly  fits  the  data,  so  making  D  go  to  one  as  pore  size 
goes  to  zero  is  best. 

Combining  the  results  from  Equations  20  and  21  into  Equation  19 
gives: 


-‘•(ft 


1 


(1+  P2y  +  P3, 


1  - 


AP 


1  +  U-  D^J 


(22) 
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Pore  size  (micro  meters) 

Figure  33  Relationship  between  parameter  D  and  pore  size 


where 


E  =  Elastic  modulus  (Young’s  modulus) 

Eo  =  Elastic  modulus  without  any  porosity 

A  =  a  constant  parameter  from  Equation  20 

L  =  pore  size  (length,  diameter,  etc.)  in  micrometers 

PI  =  0.001865 

P2  =  -7.51 

P3  =  0.8915 

P  =  volume  porosity 

Equation  22  is  good  for  aluminum  nitride  with  pore  sizes  between  5 
and  15  micrometers  since  that’s  the  data  that  was  used  to  generate  it. 
However,  to  be  useful  it  must  apply  to  other  materials  and  pore  sizes. 
Applying  Equation  22  to  cobalt  oxide’s  8  micrometer  pores  (and  letting  Eo 
=  26.1)  gives  the  results  shown  in  Figure  34,  these  results  match  the 
data  but  not  very  closely.  Figure  35  shows  the  results  of  Equation  22 
when  applied  to  silicon  carbide  (2  micrometer  pores).  These  results  are 
quite  bad  and  to  get  them  even  that  close  one  must  use  an  Eo  of  400 
Giga  Pascals,  not  the  actual  440  Giga  Pascals.. 

Linear  Derived  Equation 

Both  Figure  34  and  35  show  that  the  derived  approximate  equation 
using  the  semi-empirical  curve  fit  didn’t  work  to  well.  Better  results 
might  be  achieved  by  using  a  linear  approximation,  because  it  is  more 
generally  applicable  even  if  it  isn’t  as  accurate  for  any  one  case. 
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0.00  5.00  10.00  15.00  20.00  25.00 

Porosity  (volumetric  percentage) 

Figure  34  Cobalt  Oxide  Elastic  Modulus  vs  Porosity 
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Using  the  linear  Equation  1,  repeated  below,  means  that  only  the 
parameter  for  the  slope  (a)  needs  to  be  modified  for  the  various  pore 
sizes. 

E-Ee(l-aP)  (1) 

where 

E  =  Elastic  modulus  (Young’s  modulus) 

Eo  =  Elastic  modulus  without  any  porosity 
a  =  slope,  an  empirical  coi  stant 
P  =  volume  fraction  porosity 

Table  3  shows  all  the  data  correlating  slopes  to  pore  size,  for  the 
cases  studied  in  this  naoer. 

Material 

Si  C,  actual  test  data 
Al  N,  actual  test  data 
A1  N,  computer  data 
Co  0,  actual  test  data 
Al  N,  computer  data 
Al  N,  computer  data 
Al  N,  actual  test  data 
Al  N,  computer  data 

Table  3  Slopes  for  various  pore  sizes 

CURVEFIT  (3)  was  used  again  to  find  an  equation  that  describes  how 
slope,  a,  changes  with  oore  size.  The  result  is  again  the  Cauchy 
distribution  equation,  this  time  with  very  different  parameters: 


Pore  Size 


Micro  m. 

Slope 

2 

2.103 

5 

2.157 

5 

2.247 

8 

2.720 

8 

2.201 

10 

2.048 

10 

1.908 

15 

1.351 

where 

a  =  parameter,  slope,  for  Equation  1 
L  =  pore  size  (length,  diameter,  etc.)  in  micrometers 
PI  =  a  constant,  in  this  case  0.003866 
P2  =  a  constant,  in  this  case  -5.988 
P3  =  a  constant,  in  this  case  0.4275 

Figure  36  shows  how  Foliation  23  is  only  an  approximation  to  the 
actual  data  given  in  Tahle  3.  Putting  the  resuli  from  Equation  23  into 
Equation  1  gives: 


i 


n  1  r  t  ,  no\2 

A  A  ^  A*  *  A  **  J 


where 

E  =  Elastic  modulus  (Young’s  modulus) 

Eo  =  Elastic,  modulus  without  anv  porosity 

L  =  pore  size  (length,  diameter,  etc.)  in  micrometers 

PI  =  a  constant,  in  this  case  0.003866 

P2  =  a  constant,  in  this  case  -5.988 

P3  =  a  constant,  in  this  csrs  0.49,75 

P  =  volume  fraction  norosit.v 
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Figure  36  Relationship  between  parameter  a  and  pore  size 


Applying  Fauatinn  24  tn  thA  ra cab  that.  have  heen  looked  at  in  this 
paper  (i.e.  cobalt  otiHa.  silicon  carbide.  and  aluminum  nitride)  give  very 
accurate  results,  shown  in  Piourps  37.  38.  39.  and  40.  While  Equation  24, 
which  estimates  elastic  modulus  based  on  pore  size  and  porosity  is  not 
perfect,  it  is  close  enough  to  give  a  material  designer  a  good  starting 
estimate  of  how  a  new  material  would  behave. 

To  verify  Eauation  24.  it  should  be  used  for  some  test  case. 
Information  on  the  pore  size  of  materials  is  scarce,  so  aluminum  nitride 
will  be  used  again,  but  this  time  for  different  pore  sizes  not  previously 
examined.  Referring  hack  to  Table  1,  data  point  2  has  1.7%  porosity,  2 
micrometer  pores  (average),  and  the  same  zero  porosity  elastic  modulus 
for  aluminum  nitride  is  about  46.7  million  psia.  Putting  this  data  into 
Equation  24  gives  an  elastic  modulus  of  45.07  million  psia,  where  the 
actual  tested  numher  is  44.5  million  osia.  This  is  only  a  1.28  percent 
error,  extremely  good  for  an  empirical  first  approximation.  Several  other 
data  points  from  Tahle  1  were  tested,  and  Table  4  below  shows  the 
results.  In  all  cases  the  error  was  was  small,  but  usually  not  ignorable. 


Sample 
(Table  1) 

Porosity 

(%) 

Pore  size 
micro  ra 

E  (tested) 
million  psia 

E  (calculated) 
million  psia 

Error 

(%) 

2 

1.7 

2.0 

44.5 

45.07 

1.28 

5 

5.4 

4.0 

38.6 

41.00 

6.22 

7 

6.0 

2.5 

37.9 

40.79 

7.62 

8 

10.5 

3.0 

34.9 

36.08 

3.38 

10 

13.2 

12.0 

36.5 

35.80 

1.92 

11 

14.0 

12.0 

36.4 

35.1 

3.57 

Table  4  Equation  24  Verification  for  Aluminum  Nitride  (AIN) 
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0.00  5.00  10.00  15.00  20.00  25.00 

Porosity  (volumetric  percentage) 

Figure  37  Cobalt  Oxide  Elastic  Modulus  vs  Porosity 
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0.00  5.00  10.00  15.00  20.00  25.00 

Porosity  (volumetric  percentage) 

Figure  38  Silicon  Carbide  Elastic  Modulus  vs  Porosity 
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0.00  5.00  10.00  15.00  20.00  25.00 

Porosity  (volumetric  percentage) 

Figure  39  Aluminum  Nitride  Elastic  Modulus  vs  Porosity 
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0.00  5.00  10.00  15.00  20.00  25.00 

Porosity  (volumetric  percentage) 

Figure  40  Aluminum  Nitride  Elastic  Modulus  vs  Porosity 


Plasticity 


The  computer  model  used  did  not  t«We  into  account  any  plasticity. 

It  onlv  used  maomsconic  linear  elastic  theory,  and  yet  gave  good  results 
even  though  a  microscopic  phenomena  was  being  modeled.  No  microscopic 
phenomena,  grain  size,  plastic  deformations,  etc.,  was  considered.  The 
accuracy  is  partially  due  to  the  fact  that  at  room  temperature,  all  three 
materials  (cobalt  oxide,  silicon  carbide,  and  aluminum  nitride)  actually  are 
linearly  elastic  and  break  before  yielding.  Figure  41,  using  data  from 
Larson  and  Adams  (13:152).  shows  this  fact  on  a  macroscopic  stress  strain 
curve  for  silicon  carbide.  It  also  shows  that  temperatures  must  rise 
above  1000  °C  for  the  elastic  modulus  to  change  significantly  from  it’s 
room  temperature  value. 

I  was  unable  to  find  any  data  on  porosity  (pore  size  and 
distribution)  for  a  plastic  or  visco-elastic  material.  Future  work  using 
this  paper’s  computer  model  on  such  a  material  is  encouraged. 
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Figure  41  Stress-strain  behavior  of  Kyocera  SC-201  SiC  (Ref  13) 


VI.  Conclusion 


Up  until  now  a  material  designer  had  to  actually  make  and  test 
several  samples  of  a  material,  and  curve  fit  the  results  to  one  of  several 
empirical  formulas  (linear,  exponential  empirical,  semi-empirical,  or  the 
power  empirical  equation)  to  see  how  a  material  behaved  as  porosity 
changed.  This  thesis  has  presented  an  alternative  to  actual  production 
and  testing. 

Materials  that  are  linearly  elastic,  like  most  ceramics,  with  average 
pore  sizes  between  about  1  and  20  micrometers  could  be  evaluated  using 
Equation  24.  This  equation  gives  valuable  information  showing  how  pore 
size  and  porosity  effect  the  materials  elastic  modulus.  For  example;  a 
designer  who  wants  to  change  the  pore  size  of  a  ceramic  material  by 
using  a  cheaper  manufacturing  technique  can  get  a  good  idea  of  what  the 
new  materials  strength  will  be  using  Equation  24. 

A  closer  prediction  of  actual  material  behavior  can  be  obtained 
through  the  application  of  the  finite  element  model  (Appendix  A).  This 
computer  model  gave  accurate  results,  with  only  slight  error  when 
compared  to  lab  test  results,  for  all  materials  analyzed.  This  finite 
element  model  can  also  be  ai  ~>lied  to  materials  that  are  not  linear  elastic 
using  other  MSC/NASTRAN  solution  algorithms,  if  the  stress-strain  curve 
for  the  zeio  porosity  material  is  known. 

The  computer  model  was  also  useful  is  showing  stress  concentrations 


about  non-circular  pores.  Especially  interesting  was  the  effects  on 
stress  near  areas  where  several  pores  meet  and  join.  As  the  number  of 
touching  pores  rose,  both  the  stress  and  the  stress  gradient  on  the 
material  near  the  pore  rose  dramatically,  but  the  stress  even  one  pore 
diameter  away  changed  very  slowly,  if  at  all. 

The  computer  model  Has  clearly  shown  the  relationship  between  pore 
size  and  how  this  effected  elastic  moduli  as  the  porosity  increased  (i.e. 
the  slope  of  the  elastic  modulus  versus  porosity  line).  As  the  pore  size 
became  larger,  the  slope  decreased.  The  slope  also  decreased  for  very 
small  pores,  about  one  micrometer  small.  This  information  was  used  to 
generate  Equation  24. 

The  finite  element  model  used  in  this  thesis  represented  a  very  small 
portion  of  material.  If  a  designer  wanted  to  model  a  real  porous 
structure  they  could  take  this  model,  define  it  as  one  primary 
superelement,  and  then  using  imaging  superelements,  make  copies  of  this 
prim.'  y  superelement  and  put  them  together  like  building  blocks  to  make 
a  much  larger  structure. 

*  . 

In  summary,  the  strength  of  a  porous  material  depends  primarily  on 
volumetric  porosity,  and  pore  size.  This  thesis  has  developed  an 
equation  and  a  computer  model  which  can  be  used  to  predict  a  material’s 
elastic  modulus  using  only  pore  size,  porosity,  and  zero  porosity  elastic 


modulus. 


Appendix  A:  Sample  MSC/NASTRAN  Input  Data  File 
for  the  Computer  Model  (PORE2-9) 


ASSIGN  0UTPUT2=P0RE2-9.0P2,  STATUS=NEW,  UNIT=12 

SOL  101  $  SUPERELEMENT  STATICS 

TIME  50 

DIAG  64 

CEND 

TITLE  =  POROUS  CUBE  OF  COBOLT  OXIDE 

LOAD  =  10 

SPC  =  20 

SEALL  =  ALL 

DISP  =  ALL 

STRESS  =  ALL 

BEGIN  BULK 

$  FIRST  THE  RESIDUAL  POINTS 


GRID, 

Ill, 

9 

0.0, 

0.0,  0.0 

GRID, 

112, 

9 

0.0, 

0.0,  20.0 

GRID, 

113, 

9 

0.0, 

0.0,  40.0 

GRID, 

114, 

9 

0.0, 

0.0,  60.0 

GRID, 

115, 

9 

0.0, 

0.0,  80.0 

GRID, 

116, 

9 

0.0, 

0.0,  100.0 

GRID, 

121, 

9 

0.0, 

20.0,  0.0 

GRID, 

122, 

9 

0.0, 

20.0.  20.0 

GRID, 

123, 

9 

0.0, 

20.0,  40.0 

GRID, 

124, 

9 

0.0. 

20.0,  60.0 

GRID, 

125, 

9 

0.0, 

20.0,  80.0 

GRID, 

126, 

« 

0.0. 

20.0.  100.0 

GRID, 

131, 

9 

0.0, 

40.0,  0.0 

GRID, 

132, 

• 

0.0. 

40.0.  20.0 

GRID, 

133, 

9 

0.0, 

40.0.  40.0 

GRID. 

134. 

• 

0.0. 

40.0.  60.0 

GRID, 

135, 

9 

0.0. 

40.0.  80. 0 

GRID. 

136. 

• 

0.0. 

40.0.  100.0 

GRID, 

141, 

9 

0.0, 

60.0.  0.0 

GRID, 

142. 

• 

0.0. 

60.0.  20.0 

GRID, 

143, 

0.0, 

60.0.  40.0 

GRID, 

144. 

• 

0.0. 

60.0.  60.0 

GRID, 

145, 

9 

0.0, 

60.0.  80.0 

GRID, 

146. 

• 

0.0. 

60.0.  100.0 

GRID, 

151, 

9 

0.0, 

80.0.  0.0 

GRID, 

152. 

• 

0.0. 

80.0.  20.0 

GRID, 

153, 

9 

0.0. 

80.0.  40.0 

GRID, 

154. 

. 

0.0. 

80.0.  60  0 

GRID, 

155, 

* 

0.0. 

80.0.  80.0 

GRID, 

156. 

0.0. 

80.0.  100.0 

GRID, 

211, 

* 

20.0 

.  0.0.  0.0 

GRID. 

212. 

• 

20.0 

0  0  90  0 

GRID, 

213, 

9 

20.0 

.  0.0.  40.0 

GRID. 

214. 

. 

20.0 

.  0.0.  60.0 

GRID. 

215. 

. 

20.0. 

o.n.  ftn.n 

GRID, 

216, 

20.0. 

0.0.  100.0 

GRID. 

221. 

• 

20.0. 

20.0.  0.0 

GRID, 

222, 

4 

20.0. 

20.0.  20.0 

GRID. 

223. 

20.0. 

20.0.  40.0 

GRID, 

224, 

1 

20.0, 

20.0.  60.0 

GRID. 

225. 

20.0. 

20.0.  80.0 

GRID, 

226, 

4 

20.0. 

20.0.  100.0 

GRID, 

231. 

20.0. 

40.0.  0.0 

GRID, 

232, 

4 

20.0, 

40.0.  20.0 

GRID. 

233. 

20.0. 

40.0.  40.0 

GRID, 

234, 

4 

20.0, 

40.0.  60.0 

GRID. 

235. 

20.0. 

40.0.  80.0 

GRID, 

236, 

t 

20.0, 

40.0.  100.0 

GRID. 

241. 

« 

20.0. 

60.0.  0.0 

GRID, 

242,. 

4 

20.0. 

60.0.  20.0 

GRID, 

243. 

20.0. 

60.0.  40.0 

GRID, 

244, 

4 

20.0, 

60.0.  60.0 

GRID, 

245. 

• 

20.0. 

60.0.  80.0 

GRID, 

246, 

4 

20.0. 

60.0.  100.0 

GRID, 

251. 

• 

20.0. 

RO.O.  0.0 

GRID, 

252, 

4 

20.0, 

80.0.  20.0 

GRID, 

253. 

« 

20.0. 

80.0.  40.0 

GRID, 

254, 

4 

20.0, 

80.0,  60.0 

GRID, 

255, 

• 

20.0. 

80.0.  80.0 

GRID, 

256, 

4 

20.0, 

80.0,  100.0 

GRID, 

311, 

4 

40.0. 

0.0,  0.0 

GRID, 

312, 

4 

40.0, 

0.0,  20.0 

GRID, 

313, 

4 

40.0, 

0.0,  40.0 

GRID, 

314, 

1 

40.0, 

0.0,  60.0 

GRID, 

315, 

1 

40.0. 

0.0.  80.0 

GRID, 

316, 

4 

40.0, 

0.0,  100.0 

GRID, 

321, 

4 

40.0, 

20.0,  0.0 

GRID, 

322, 

4 

40.0, 

20.0,  20.0 

GRID, 

323, 

4 

40.0, 

20.0,  40.0 

GRID, 

324, 

4 

40.0, 

20.0,  60.0 

GRID, 

325, 

4 

40.0, 

20.0,  80.0 

GRID, 

326, 

4 

40.0, 

20.0,  100.0 

GRID, 

331, 

4 

40.0, 

40.0,  0.0 

GRID, 

332, 

4 

40.0, 

40.0,  20.0 

GRID, 

333, 

4 

40.0, 

40.0,  40.0 

GRID, 

334, 

4 

40.0, 

40.0,  60.0 

GRID, 

335, 

4 

40.0, 

40.0,  80.0 

GRID, 

336, 

4 

40.0, 

40.0,  100.0 

GRID, 

341, 

4 

40.0, 

60.0,  0.0 

GRID, 

342, 

4 

40.0, 

60.0,  20.0 

GRID, 

343, 

4 

40.0, 

60.0,  40.0 

GRID, 

344, 

4 

40.0, 

60.0,  60.0 

GRID, 

345, 

4 

40.0, 

60.0,  80.0 

GRID, 

346, 

4 

40.0, 

60.0,  100.0 

GRID, 

351, 

4 

40.0, 

80.0,  0.0 

GRID,  352,  ,  40.0.  80.0.  90  o 
GRID,  353,  ,  40.0,  80.0,  40.0 
GRID,  354,  ,  40.0.  80.0,  60.0 
GRID,  355,  ,  40.0,  80.0,  80.0 
GRID,  356,  .  40.0.  80.0.  100.0 
GRID,  411,  ,  60.0,  0.0,  0.0 
GRID,  412.  .  60.0.  0.0.  20.0 
GRID,  413,  ,  60.0,  0.0.  40.0 
GRID,  414.  .  60.0.  0.0.  60.0 
GRID,  415,  ,  60.0,  0.0,  80.0 
GRID,  416,  ,  60.0.  0.0.  100.0 
GRID,  421,  ,  60.0,  20.0,  0.0 
GRID,  422,  ,  60.0,  20.0,  20.0 
GRID,  423,  ,  60.0,  20.0,  40.0 
GRID,  424,  ,  60.0,  20.0,  60.0 
GRID,  425,  ,  60.0,  20.0,  80.0 
GRID,  426,  ,  60.0,  20.0.  100.0 
GRID,  431,  ,  60.0,  40.0,  0.0 
GRID,  432,  ,  60.0,  40.0,  20.0 
GRID,  433,  ,  60.0,  40.0,  40.0 
GRID,  434,  ,  60.0.  40.0.  60.0 
GRID,  435,  ,  60.0,  40.0,  80.0 
GRID,  436,  ,  60.0.  40.0,  100.0 
GRID,  441,  ,  60.0,  60.0,  0.0 
GRID,  442,  ,  60.0.  60.0.  20.0 
GRID,  443,  ,  60.0,  60.0,  40.0 
GRID,  444,  ,  60.0.  60.0,  60.0 
GRID,  445,  ,  60.0,  60.0,  80.0 
GRID,  446,  ,  60.0.  60.0,  100.0 
GRID,  451,  ,  60.0,  80.0,  0.0 
GRID,  452.  .  60.0.  80.0.  20.0 
GRID,  453,  ,  60.0,  80.0.  40.0 
GRID,  454.  .  60.0.  80.0.  60.0 
GRID,  455,  ,  60.0,  80.0.  80.0 
GRID,  456.  .  60.0.  80.0,  100.0 
GRID,  511,  ,  80.0,  0.0,  0.0 
GRID,  512,  ,  80.0.  0.0.  20.0 
GRID,  513,  ,  80.0,  0.0.  40.0 
GRID,  514.  .  80.0.  0.0.  60.0 
GRID,  515,  ,  80.0,  0.0,  80.0 
GRID,  516,  .  80.0.  0.0.  100.0 
GRID,  521,  ,  80.0,  20.0,  0.0 
GRID,  522.  .  80.0.  20.0.  20.0 
GRID,  523,  ,  80.0,  20.0,  40.0 
GRID,  524,  ,  80.0,  20.0,  60.0 
GRID,  525,  ,  80.0,  20.0,  80.0 
GRID,  526,  ,  80.0,  20.0,  100.0 
GRID,  531,  ,  80.0,  40.0,  0.0 
GRID,  532,  ,  80.0,  40.0,  20.0 
GRID,  533,  ,  80.0,  40.0,  40.0 
GRID,  534,  ,  80.0,  40.0,  60.0 
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GRID,  535,  ,  80.0,  40.0,  80.0 
GRID,  536,  ,  80.0,  40.0,  100.0 
GRID,  541,  ,  80.0,  60.0,  0.0 
GRID,  542,  ,  80.0,  60.0,  20.0 
GRID,  543,  ,  80.0,  60.0,  40.0 
GRID,  544,  ,  80.0,  60.0,  60.0 
GRID,  545,  ,  80.0,  60.0,  80.0 
GRID,  546,  ,  80.0,  60.0,  100.0 
GRID,  551,  ,  80.0,  80.0,  0.0 
GRID,  552,  ,  80.0,  80.0,  20.0 
GRID,  553,  ,  80.0,  80.0,  40.0 
GRID,  554,  ,  80.0,  80.0,  60.0 
GRID,  555,  ,  80.0,  80.0,  80.0 
GRID,  556,  ,  80.0,  80.0,  100.0 
$  ************************************************** 
$  *  SINGLE  CELL,  SOLID 
GRID,  1011,  ,  66.0,  6.0,  26.0 
GRID,  1012,  ,  74.0,  6.0,  26.0 
GRID,  1013,  ,  74.0,  6.0,  34.0 
GRID,  1014,  ,  66.0,  6.0,  34.0 
GRID,  1015,  ,  74.0,  14.0,  26.0 
GRID,  1016,  ,  74.0,  14.0,  34.0 
GRID,  1017,  ,  66.0,  14.0,  34.0 
GRID,  1018,  ,  66.0,  14.0,  26.0 
SESET,  20,  1011,  THRU,  1018 
CHEXA,  1,  30,  412,  512,  513,  413,  1011,  1012 
,  1013,  1014 

CHEXA,  2,  30,  512,  522,  523,  513,  1012,  1015 
,  1016,  1013 

CHEXA,  3,  30,  522,  523,  423,  422,  1015,  1016 
,  1017,  1018 

CHEXA,  4,  30,  422,  412,  413,  423,  1018.  1011 
,  1014,  1017 

CHEXA,  5,  30,  413,  513.  523.  423.  1014.  1013 
,  1016,  1017 

CHEXA,  6,  30,  412,  512.  522.  422.  1011.  1012 
,  1015.  1018 

CHEXA,  7,  30,  1011,  1012.  1013.  1014.  1018.  1015 
,  1016.  1017 
PSOLID,  30,  40,  .  3 
MAT1,  40.  2.0202+11.  .  0.2Q5.  1.0 
$  ************************************************** 
$  *  SINGLE  CELL.  PORE 
GRID,  1001,  .  6.0,  6.0,  6.0 
GRID,  1002.  .  14.0.  fi.O.  fi.O 
GRID,  1003,  ,  14.0,  6.0.  14.0 
GRID,  1004.  .  6.0.  6.0.  14.0 
GRID,  1005,  ,  14.0,  14.0.  6.0 
GRID,  1006.  .  14.0.  14.0.  14.0 
GRID,  1007,  ,  6.0,  14.0.  14.0 
GRID,  1008.  .  6.0.  14.0.  6.0 


SESET.  1.  1001.  THRU.  1008 

CHEXA,  11,  30,  111,  211,  212.  112.  1001.  1002 

,  1003,  1004 

CHEXA,  12,  30,  211,  221,  222.  212.  1002.  1005 
,  1006,  1003 

CHEXA,  13,  30,  221,  222,  122,  121.  1005.  1006 
,  1007,  1008 

CHEXA,  14,  30,  121,  111.  112.  122.  1008.  1001 
,  1004,  1007 

CHEXA,  15,  30,  112,  212,  222,  122.  1004.  1003 
,  1006,  1007 

CHEXA,  16,  30,  111,  211.  221.  121,  1001,  1002 
,  1005,  1008 

CHEXA,  17,  31,  1001,  1002,  1003,  1004,  1008,  1005 
,  1006,  1007 

$  ************************************************** 

$  *  SINGLE  CELL,  TWO  PORES 

GRID,  2001,  ,  61.5,  6.0,  6.0 

GRID,  2002,  ,  69.5,  6.0,  6.0 

GRID,  2003,  ,  69.5,  6.0,  14.0 

GRID,  2004,  ,  61.5,  6.0,  14.0 

GRID,  2005,  ,  69.5,  14.0,  6.0 

GRID,  2006,  ,  69.5,  14.0,  14.0 

GRID,  2007,  ,  61.5,  14.0,  14.0 

GRID,  2008,  ,  61.5,  14.0,  6.0 

GRID,  2009,  ,  70.5,  6.0,  6.0 

GRID,  2010,  ,  78.5,  6.0,  6.0 

GRID,  2011,  ,  78.5,  6.0,  14.0 

GRID,  2012,  ,  70.5,  6.0,  14.0 

GRID,  2013,  ,  78.5,  14.0,  6.0 

GRID,  2014,  ,  78.5.  14.0.  14.0 

GRID,  2015,  ,  70.5,  14.0,  14.0 

GRID,  2016,  ,  70.5.  14.0.  6.0 

SESET,  4,  2001,  THRU,  2016 

CHEXA,  20,  30.  411.  511.  512.  412.  2002.  2009 

,  2012,  2003 

CHEXA,  21,  30,  511,  521.  522.  512.  2010.  2013 
,  2014,  2011 

CHEXA,  22,  30,  521.  421.  422.  522.  2016.  2005 
,  2006,  2015 

CHEXA,  23,  30,  421,  411.  412.  422.  2008.  2001 
,  2004,  2007 

CHEXA,  24,  30,  412,  512.  522.  422.  2003.  2012 
,  2015,  2006 

CHEXA,  25,  30,  411,  511.  521.  421.  2002.  2009 
,  2016,  2005 

CHEXA,  26,  30.  2002.  2009.  2012.  2003.  2005.  2016 
,  2015.  2006 

CHEXA,  27,  31,  2001.  2002.  2003.  2004.  2008.  2005 
,  2006.  2007 

CHEXA,  28,  31,  2009.  2010.  2011.  2012.  2016.  2013 


,  2014.  2015 

CPENTA,  29.  30.  411.  2001.  2002.  412.  9.004.  2003 

CPENTA.  30.  30.  511.  2010.  200ft.  519.  9.011.  2012 

CPENTA,  31,  30.  521.  2013.  201  fi.  522.  2014.  2015 

CPENTA.  32.  30.  421.  200R.  2005.  422.  9007.  200fi 

CPENTA,  33.  30.  512.  2011.  9.012.  522.  2014.  2015 

CPENTA.  34.  30.  412.  2004.  2003.  422.  2007.  200fi 

CPENTA,  35,  30.  411.  2001.  2002.  421.  2.00ft.  2005 

CPENTA.  3fi.  30.  511.  2010.  200ft.  521.  2013.  2015 

$  ******************************************** 

$  SINGLE  CELL.  THREE  PORES 

GRID,  3111,  ,  2.0.  2.0.  42.0 

GRID,  3112.  .  2.0.  2.0.  50.0 

GRID,  3113,  .  2.0.  2.0.  5fi.O 

GRID,  3121.  .  2.0.  10.0.  42.0 

GRID,  3122,  ,  2.0.  10.0.  50.0 

GRID.  3123.  .  2.0.  10.0.  5ft.O 

GRID,  3131,  .  2.0.  Ift.O.  42.0 

GRID,  3132.  .  2.0.  Ift.O.  50.0 

GRID,  3133,  .  2.0.  18.0.  58.0 

GRID,  3211.  .  10.0.  2.0.  42.0 

GRID,  3212,  ,  10.0.  2.0.  50.0 

GRID.  3213.  .  10.0.  2.0.  5R.0 

GRID,  3221,  ,  10.0,  10.0.  42.0 

GRID,  3222.  .  10.0.  10.0.  50.0 

GRID,  3223,  ,  10.0.  10.0.  58.0 

GRID,  3231.  .  10.0.  Ift.O.  49.0 

GRID,  3232,  ,  10.0.  18.0.  50.0 

GRID.  3233.  .  10.0.  Ift.O.  5R.0 

GRID,  3311,  ,  18.0.  2.0.  42.0 

GRID.  3312.  .  18.0.  2.0.  50.0 

GRID,  3313,  ,  18.0.  2.0.  58.0 

GRID,  3321.  .  18.0.  10.0.  42.0 

GRID,  3322,  ,  18.0,  10.0.  50.0 

GRID,  3323.  .  18.0.  10.0.  58.0 

GRID,  3331,  ,  18.0,  18.0,  42.0 

GRID,  3332,  ,  18.0,  18.0,  50.0 

GRID,  3333,  ,  18.0,  18.0,  58.0 

CHEXA,  40,  30,  113,  213,  214,  114,  3111,  3311 

,  3313,  3113 

CHEXA,  41,  30,  213,  223,  224,  214,  3311,  3331 
,  3333,  3313 

CHEXA,  42,  30,  223,  123,  124,  224,  3331,  3131 
,  3133,  3333 

CHEXA,  43,  30,  123,  113,  114,  124,  3131,  3111 
,  3113,  3133 

CHEXA,  44,  30,  114,  214,  224,  124,  3113,  3313 
,  3333,  3133 

CHEXA,  45,  30,  113,  213,  223,  123,  3111,  3311 
,  3331,  3131 

CHEXA,  46,  31,  3111,  3211.  3221.  3121.  3112.  3212 


,  3222,  3122 

CHEXA,  47,  30,  3211,  3311.  3321.  3221.  3212.  3312 
,  3322,  3222 

CHEXA,  48,  30,  3121,  3221.  3231.  3131.  3122.  3222 
,  3232.  3132 

CHEXA,  49,  30,  3221,  3321.  3331.  32.31.  3222.  3322 
,  3332,  3232 

CHEXA,  50,  30,  3112.  3212.  3222.  3122.  3113.  3213 
,  3223.  3123 

CHEXA,  51,  31,  3212,  3312,  3322.  3222.  3213.  3313 
,  3323.  3223 

CHEXA,  52,  31,  3122.  3222.  .32.32.  3132.  3123.  3223 
,  3233,  3133 

CHEXA,  53,  30,  3222.  3322.  .3.332.  3232.  .322.3.  3.323 
,  3333,  3233 

SESET,  33,  3111,  THRU,  3333 

$  *:M:***************************************** 

$  SINGLE  CELL,  FOUR  PORES 

GRID,  4111.  .  2.0.  2.0.  62.0 

GRID,  4112,  ,  2.0,  2.0,  70.0 

GRID,  4113,  .  2.0.  2.0.  78.0 

GRID,  4121,  ,  2.0,  10.0,  62.0 

GRID,  4122,  ,  2.0,  10.0,  70.0 

GRID,  4123,  ,  2.0,  10.0,  78.0 

GRID,  4131,  ,  2.0,  18.0,  62.0 

GRID,  4132,  ,  2.0,  18.0,  70.0 

GRID,  4133,  ,  2.0,  18.0,  78.0 

GRID,  4211,  ,  10.0,  2.0,  62.0 

GRID,  4212,  ,  10.0,  2.0,  70.0 

GRID,  4213,  ,  10.0,  2.0,  78.0 

GRID,  4221,  ,  10.0,  10.0,  62.0 

GRID,  4222,  ,  10.0,  10.0,  70.0 

GRID,  4223,  ,  10.0,  10.0,  78.0 

GRID,  4231,  ,  10.0,  18.0,  62.0 

GRID,  4232,  ,  10.0,  18.0,  70.0 

GRID,  4233,  ,  10.0,  18.0,  78.0 

GRID,  4311,  ,  18.0,  2.0,  62.0 

GRID,  4312,  ,  18.0,  2.0,  70.0 

GRID,  4313,  ,  18.0,  2.0,  78.0 

GRID,  4321,  ,  18.0,  10.0,  62.0 

GRID,  4322,  ,  18.0,  10.0,  70.0 

GRID,  4323,  ,  18.0,  10.0,  78.0 

GRID,  4331,  ,  18.0,  18.0,  62.0 

GRID,  4332,  ,  18.0,  18.0,  70.0 

GRID,  4333,  ,  18.0,  18.0,  78.0 

CHEXA,  60,  30,  114,  214,  215,  115,  4111,  4311 

,  4313,  4113 

CHEXA,  61,  30,  214,  224,  225,  215.  4311.  4331 
,  4333,  4313 

CHEXA,  62,  30,  224,  124,  125,  225,  4331,  4131 
,  4133,  4333 


CHEXA,  63,  30,  124,  114,  115,  125,  4131,  4111 
,  4113,  4133 

CHEXA,  64,  30,  115,  215,  225,  125,  4113,  4313 
4333  4133 

CHEXA,  65,  30,  114,  214,  224,  124,  4111,  4311 
,  4331,  4131 

CHEXA,  66,  31,  4111,  4211,  4221,  4121,  4112,  4212 
,  4222,  4122 

CHEXA,  67,  30,  4211,  4311,  4321,  4221,  4212.  4312 
4322  4222 

CHEXa’,  68,  30,  4121,  4221.  4231.  4131.  4122.  4222 
,  4232,  4132 

CHEXA,  69,  31,  4221.  4321.  4331.  4231.  4222.  4322 
4332  4232 

CHEXA,  54,  30,  4112,  4212.  4222.  4122.  4113.  4213 
,  4223.  4123 

CHEXA,  55,  31,  4212,  4312.  4322.  4222.  4213.  4313 
,  4323,  4223 

CHEXA,  56,  31,  4122,  4222,  4232.  4132.  4123.  4223 
,  4233,  4133 

CHEXA,  57,  30,  4222,  4322,  4332,  4232,  4223,  4323 
,  4333,  4233 

SESET,  49,  4111,  THRU,  4333 
$  ******************************************** 

$  SINGLE  CELL,  FIVE  PORES 


GRID, 

5111, 

9 

2.0,  62.0,  42.0 

GRID, 

5112, 

9 

2.0,  62.0,  50.0 

GRID, 

5113, 

9 

2.0,  62.0,  58.0 

GRID, 

5121, 

9 

2.0,  70.0,  42.0 

GRID, 

5122, 

9 

2.0,  70.0,  50.0 

GRID, 

5123, 

9 

2.0,  70.0,  58.0 

GRID, 

5131, 

9 

2.0,  78.0,  42.0 

GRID, 

5132, 

9 

2.0,  78.0,  50.0 

GRID, 

5133, 

9 

2.0,  78.0,  58.0 

GRID, 

5211, 

9 

10.0,  62.0,  42.0 

GRID, 

5212, 

9 

10.0,  62.0,  50.0 

GRID, 

5213, 

9 

10.0,  62.0,  58.0 

GRID, 

5221, 

9 

10.0,  70.0,  42.0 

GRID, 

5222, 

9 

10.0,  70.0,  50.0 

GRID, 

5223, 

9 

10.0,  70.0,  58.0 

GRID, 

5231, 

9 

10.0,  78.0,  42.0 

GRID, 

5232, 

9 

10.0,  78.0,  50.0 

GRID, 

5233, 

9 

10.0,  78.0,  58.0 

GRID, 

5311, 

9 

18.0,  62.0,  42.0 

GRID, 

5312, 

9 

18.0,  62.0,  50.0 

GRID, 

5313, 

9 

18.0,  62.0,  58.0 

GRID, 

5321, 

9 

18.0,  70.0,  42.0 

GRID, 

5322, 

9 

18.0,  70.0,  50.0 

GRID, 

5323, 

9 

18.0,  70.0,  58.0 

GRID, 

5331, 

9 

18.0,  78.0,  42.0 

GRID, 

5332, 

9 

18.0,  78.0,  50.0 
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GRID,  5333,  ,  18.0,  78.0,  58.0 

CHEXA,  70,  30,  143,  243,  244,  144,  5111,  5311 

,  5313,  5113 

CHEXA,  71,  30,  243,  253,  254,  244,  5311,  5331 
,  5333,  5313 

CHEXA,  72,  30,  253,  153,  154,  254,  5331.  5131 
5133  5333 

CHEXA,  73,  30,  153,  143,  144,  154,  5131,  5111 
,  5113,  5133 

CHEXA,  74,  30,  144,  244,  254,  154,  5113,  5313 
,  5333,  5133 

CHEXA,  75,  30,  143,  243,  253,  153,  5111,  5311 
,  5331,  5131 


• 

CHEXA,  76,  31, 

,  5222,  5122 

5111, 

5211, 

5221, 

5121, 

5112, 

5212 

CHEXA,  77,  31, 

,  5322,  5222 

5211, 

5311, 

5321, 

5221, 

5212, 

5312 

CHEXA,  78,  30, 

,  5232,  5132 

5121, 

5221, 

5231, 

5131, 

5122, 

5222 

• 

CHEXA,  79,  31, 

,  5332,  5232 

5221, 

5321, 

5331, 

5231, 

5222, 

5322 

CHEXA,  80,  30, 

,  5223,  5123 

5112, 

5212, 

5222, 

5122, 

5113, 

5213 

CHEXA,  81,  31, 

,  5323,  5223 

5212, 

5312, 

5322, 

5222, 

5213, 

5313 

• 

CHEXA,  82,  31, 

,  5233,  5133 

5122, 

5222, 

5232, 

5132, 

5123, 

5223 

CHEXA,  83,  30, 

5222, 

5322, 

5332, 

5232, 

5223, 

5323 

,  5333,  5233 

SESET,  45,  5111,  THRU,  5333 

$  ;Mt************:M:**************************** 

$  SINGLE  CELL,  SIX  PORES 


GRID, 

6111, 

9 

62.0, 

2.0,  42.0 

GRID, 

6112, 

• 

62.0, 

2.0,  50.0 

GRID, 

6113, 

9 

62.0, 

2.0,  58.0 

GRID, 

6121, 

9 

62.0, 

10.0,  42.0 

GRID, 

6122, 

1 

62.0, 

10.0,  50.0 

GRID, 

6123, 

9 

62.0, 

10.0,  58.0 

GRID, 

6131, 

9 

62.0, 

18.0,  42.0 

GRID, 

6132, 

9 

62.0, 

18.0,  50.0 

GRID, 

6133, 

9 

62.0, 

18.0,  58.0 

GRID, 

6211, 

9 

70.0, 

2.0,  42.0 

GRID, 

6212, 

9 

70.0, 

2.0,  50.0 

GRID, 

6213, 

9 

70.0, 

2.0,  58.0 

GRID, 

6221, 

9 

70.0, 

10.0,  42.0 

GRID, 

6222, 

9 

70.0, 

10.0,  50.0 

GRID, 

6223, 

9 

70.0, 

10.0,  58.0 

GRID, 

6231, 

9 

70.0, 

18.0,  42.0 

GRID, 

6232, 

9 

70.0, 

18.0,  50.0 

GRID, 

6233, 

9 

70.0, 

18.0,  58.0 

GRID, 

6311, 

9 

78.0, 

2.0,  42.0 
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GRID,  6312,  ,  78.0,  2.0,  50.0 

GRID,  6313,  ,  78.0,  2.0,  58.0 

GRID,  6321,  ,  78.0,  10.0,  42.0 

GRID,  6322,  ,  78.0,  10.0,  50.0 

GRID,  6323,  ,  78.0,  10.0,  58.0 

GRID,  6331,  ,  78.0,  18.0,  42.0 

GRID,  6332,  ,  78.0,  18.0,  50.0 

GRID,  6333,  ,  78.0,  18.0,  58.0 

CHEXA,  90,  30.  413.  513.  514.  414.  6111.  6311 

,  6313,  6113 

CHEXA,  91,  30,  513,  523,  524.  514.  6311.  6331 
,  6333,  6313 

CHEXA,  92,  30,  523,  423.  424.  524.  6331.  6131 
,  6133,  6333 

CHEXA,  93,  30,  423.  413.  414.  424.  6131.  6111 
,  6113,  6133 

CHEXA,  94,  30,  414,  514,  524,  424,  6113,  6313 
,  6333,  6133 

CHEXA,  95,  30,  413,  513,  523,  423,  6111,  6311 
,  6331,  6131 

CHEXA,  96,  31,  6111,  6211,  6221,  6121,  6112,  6212 
,  6222,  6122 

CHEXA,  97,  31,  6211,  6311,  6321,  6221,  6212,  6312 
,  6322,  6222 

CHEXA,  98,  30,  6121,  6221,  6231,  6131,  6122,  6222 
,  6232,  6132 

CHEXA,  99,  31,  6221,  6321,  6331,  6231,  6222,  6322 
,  6332,  6232 

CHEXA,  84,  31,  6112,  6212,  6222,  6122,  6113,  6213 
,  6223,  6123 

CHEXA,  85,  30,  6212,  6312,  6322,  6222,  6213,  6313 
,  6323,  6223 

CHEXA,  86,  31,  6122,  6222,  6232,  6132,  6123,  6223 
,  6233,  6133 

CHEXA,  87,  30,  6222,  6322,  6332,  6232,  6223,  6323 
,  6333,  6233 

SESET,  36,  6111,  THRU,  6333 
$  ************^*****:*:**********#************* 

$  SINGLE  CELL,  SEVEN  PORES 
GRID,  7111,  ,  2.0,  2.0,  22.0 
GRID,  7112,  ,  2.0,  2.0,  30.0 
GRID,  7113,  ,  2.0,  2.0,  38.0 
GRID,  7121,  ,  2.0,  10.0,  22.0 
GRID,  7122,  ,  2.0,  10.0,  30.0 
GRID,  7123,  ,  2.0,  10.0,  38.0 
GRID,  7131,  ,  2.0,  18.0,  22.0 
GRID,  7132,  ,  2.0,  18.0,  30.0 
GRID,  7133,  ,  2.0,  18.0,  38.0 
GRID,  7211,  ,  10.0,  2.0,  22.0 
GRID,  7212,  ,  10.0,  2.0,  30.0 
GRID,  7213,  ,  10.0,  2.0,  38.0 
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GRID,  7221,  , 
GRID,  7222,  , 
GRID,  7223,  , 
GRID,  7231,  , 


GRID,  7232, 


10.0,  10.0,  22.0 
10.0,  10.0,  30.0 
10.0,  10.0,  38.0 
10.0,  18.0,  22.0 
10.0,  18.0,  30.0 
10.0,  18.0,  38.0 
18.0,  2.0,  22.0 
18.0,  2.0,  30.0 
18.0,  2.0,  38.0 
18.0,  10.0,  22.0 
GRID,  7322,  ,  18.0,  10.0.  30.0 
GRID,  7323,  ,  18.0,  10.0,  38.0 
GRID,  7331,  ,  18.0,  18.0,  22.0 
GRID,  7332,  ,  18.0,  18.0,  30.0 
GRID,  7333,  ,  18.0,  18.0,  38.0 


GRID,  7233,  , 
GRID,  7311,  , 
GRID,  7312,  , 


GRID,  7313, 


GRID,  7321, 


CHEXA,  170, 

,  7313,  7113 

CO 

o 

112, 

212, 

213, 

113, 

7111, 

7311 

CHEXA,  171, 

,  7333,  7313 

30, 

212, 

222, 

223, 

213, 

7311, 

7331 

CHEXA,  172, 

,  7133,  7333 

30, 

222, 

122, 

123, 

223, 

7331, 

7131 

CHEXA,  173, 

,  7113,  7133 

30, 

122, 

112, 

113, 

123, 

7131, 

7111 

CHEXA,  174, 

,  7333,  7133 

30, 

113, 

213, 

223, 

123, 

7113, 

7313 

CHEXA,  175, 

,  7331,  7131 

o 

CO 

112, 

212, 

222, 

122, 

7111, 

7311 

CHEXA,  176,  31,  7111,  7211,  7221,  7121,  7112,  7212 
,  7222,  7122 

CHEXA,  177,  31,  7211,  7311,  7321,  7221,  7212,  7312 
,  7322,  7222 

CHEXA,  178,  31,  7121,  7221,  7231,  7131,  7122,  7222 
,  7232,  7132 

CHEXA,  179,  31,  7221,  7321,  7331,  7231,  7222.  7322 
,  7332,  7232 

CHEXA,  180,  30,  7112,  7212,  7222,  7122,  7113,  7213 
,  7223,  7123 

CHEXA,  181,  31,  7212,  7312,  7322,  7222,  7213,  7313 
,  7323,  7223 

CHEXA,  182,  31,  7122,  7222,  7232,  7132,  7123,  7223 
,  7233,  7133 

CHEXA,  183,  31,  7222,  7322,  7332,  7232,  7223,  7323 
,  7333,  7233 

SESET,  17,  7111,  THRU,  7333 
$  ******************************************** 

$  SINGLE  CELL,  EIGHT  PORE  CONGLOMERATE 

GRID,  8111,  ,  22.0,  22.0,  22.0 

GRID,  8113,  ,  22.0,  22.0,  38.0 

GRID,  8131,  ,  22.0,  38.0,  22.0 

GRID,  8133,  ,  22.0,  38.0,  38.0 

GRID,  8311,  ,  38.0,  22.0,  22.0 
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GRID,  8313,  ,  38.0,  22.0,  38.0 

GRID,  8331,  ,  38.0,  38.0,  22.0 

GRID,  8333,  ,  38.0,  38.0,  38.0 

CHEXA,  8,  30,  222,  322,  323,  223,  8111,  8311 

,  8313,  8113 

CHEXA,  9,  30,  322,  332,  333,  323,  8311,  8331 
,  8333,  8313 

CHEXA,  18,  30,  332,  232,  233,  333,  8331,  8131 
,  8133,  8333 

CHEXA,  38,  30,  232,  222,  223,  233,  8131.  8111 
,  8113,  8133 

CHEXA,  58,  30,  223,  323,  333,  233,  8113.  8313 
,  8333,  8133 

CHEXA,  88,  30,  222,  322,  332,  232,  8111.  8311 
,  8331,  8131 

CHEXA,  10,  31,  8111,  8311,  8313,  8113.  8131.  8331 
,  8333,  8133 

SESET,  38,  8111,  THRU,  8333 
$  ****  COPIED  SUPERFT.EMENTR 
CSUPER,  2,  45,  211,  212,  221,  222,  311,  312 
,  321,  322 

C SUPER,  3,  4,  311,  312,  321,  322,  411,  412 
,  421,  422 

CSUPER,  5,  49,  121,  122,  131,  132,  221,  222 
,  231,  232 

CSUPER,  6,  33,  221,  222,  231,  232,  321,  322 
,  331,  332 

CSUPER,  7,  4,  321,  322,  331,  332,  421,  422 
,  431,  432 

CSUPER,  8,  36,  421,  422,  431,  432,  521,  522 
,  531,  532 

CSUPER,  9,  33,  131,  132,  141,  142,  231,  232 
,  241,  242 

CSUPER,  10,  4,  231,  232,  241,  242,  331,  332 
,  341,  342 

CSUPER,  11,  33,  331,  332,  341,  342,  431,  432 
,  441,  442 

CSUPER,  12,  49,  431,  432,  441,  442,  531,  532 
541  542 

CSUPER,  13,  33,  141,  142,  151,  152,  241,  242 
251  252 

CSUPER,  14,  49,  241,  242,  251,  2o2,  341,  342 
,  351,  352 

CSUPER,  15,  4,  341,  342,  351,  352,  441,  442 
,  451,  452 

CSUPER,  16,  33,  441,  442,  451,  452,  541,  542 
,  551,  552 

CSUPER,  18,  4,  212,  213,  222,  223,  312,  313 
,  322,  323 

CSUPER,  19,  49,  312,  313,  322,  323,  412,  413 
,  422,  423 


CSUPER,  21,  4,  122,  123,  132,  133,  222,  223 
,  232,  233 

CSUPER,  22,  33,  222,  223,  232,  233,  322,  323 
,  332,  333 

CSUPER,  23,  49,  322,  323,  332,  333,  422,  423 
,  432,  433 

CSUPER,  24,  4,  422,  423,  432,  433,  522,  523 
,  532,  533 

CSUPER,  25,  33,  132,  133,  142,  143,  232,  233 
,  242,  243 

CSUPER,  26,  33,  232,  233,  242,  243,  332,  333 
,  342,  343 

CSUPER,  27,  1,  332,  333,  342,  343,  432,  433 
,  442,  443 

CSUPER,  28,  45,  432,  433,  442,  443,  532,  533 
,  542,  543 

CSUPER,  29,  4,  142,  143,  152,  153,  242,  243 
,  252,  253 

CSUPER,  30,  33,  242,  243,  252,  253,  342,  343 
,  352,  353 

CSUPER,  31,  45,  342,  343,  352,  353,  442,  443 
452  453 

CSUPER,  32,  49,  442,  443,  452,  453,  542,  543 
,  552,  553 

CSUPER,  34,  49,  213,  214,  223,  224,  313,  314 
323  324 

CSUPER,  35,  49,  313,  314,  323,  324,  413,  414 
,  423,  424 

CSUPER,  37,  45,  123,  124,  133,  134,  223,  224 
,  233,  234 

CSUPER,  39,  36,  323,  324,  333,  334,  423,  424 
,  433,  434 

CSUPER,  40,  4,  423,  424,  433,  434,  523,  524 
,  533,  534 

CSUPER,  41,  45,  133,  134,  143,  144,  233,  234 
,  243,  244 

CSUPER,  42,  4,  233,  234,  243,  244,  333,  334 
,  343,  344 

CSUPER,  43,  49,  333,  334,  343,  344,  433,  434 
,  443,  444 

CSUPER,  44,  4,  433,  434,  443,  444,  533,  534 
,  543,  544 

CSUPER,  46,  17,  243,  244,  253,  254.  343.  344 
,  353,  354 

CSUPER,  47,  33,  343,  344,  353,  354,  443,  444 
,  453,  454 

CSUPER,  48,  49,  443,  444,  453,  454,  543,  544 
,  553,  554 

CSUPER,  50,  45,  214,  215,  224,  225,  314,  315 
,  324,  325 

CSUPER,  51,  1,  314,  315,  324,  325,  414,  415 


,  424,  425 
CSUPER,  52, 
,  524,  525 
CSUPER,  53, 
,  234,  235 
CSUPER,  54, 
,  334,  335 
CSUPER,  55, 
,  434,  435 
CSUPER,  56, 
,  534,  535 
CSUPER,  57, 
,  244,  245 
CSUPER,  58, 
,  344,  345 
CSUPER,  59, 
,  444,  445 
CSUPER,  60, 
,  544,  545 
CSUPER,  61, 
,  254,  255 
CSUPER,  62, 
,  354,  355 
CSUPER,  63, 
,  454,  455 
CSUPER,  64, 
,  554,  555 
CSUPER,  65, 
,  225,  226 
CSUPER,  66, 
,  325,  326 
CSUPER,  67, 
,  425,  426 
CSUPER,  68, 
,  525,  526 
CSUPER,  69, 
,  235,  236 
CSUPER,  70, 
,  335,  336 
CSUPER,  71, 
,  435,  436 
CSUPER,  72, 
,  535,  5oo 
CSUPER,  73, 
,  245,  246 
CSUPER,  74, 
,  345,  346 
CSUPER,  75, 
,  445,  446 
CSUPER,  76, 
,  545,  546 


45,  414,  415,  424,  425,  514,  515 
4,  124,  125,  134,  135,  224,  225 
36,  224,  225,  234,  235,  324,  325 

33,  324,  325,  334,  335,  424,  425 

1,  424,  425,  434,  435,  524,  525 
49,  134,  135,  144,  145,  234,  235 

33,  234,  235,  244,  245,  334,  335 

33,  334,  335,  344,  345,  434,  435 

1,  434,  435,  444,  445,  534,  535 
45,  144,  145,  154,  155,  244,  245 
4,  244,  245,  254,  255,  344,  345 
33,  344,  345,  354,  355,  444,  445 
1,  444,  445,  454,  455,  544,  545 
33,  115,  116,  125,  126,  215,  216 

20,  215,  216,  225,  226,  315,  316 

49,  315,  316,  325,  326,  415,  416 

33,  415,  416,  425,  426,  515,  516 

49,  125,  126,  135,  136,  225,  226 

1,  225,  226,  235,  236,  325,  326 
33,  325,  326,  335,  336,  425,  426 
4,  425,  426,  435,  436,  525,  526 
33,  135,  136,  145,  146,  235,  236 

20,  235,  236,  245,  246,  335,  336 

4,  335,  336,  345,  346,  435,  436 
49,  435,  436,  445,  446,  535,  536 
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CSUPER,  77,  33,  145,  146,  155,  156,  245,  246 
,  255,  256 

CSUPER,  78,  1,  245,  246,  255,  256,  345,  346 
,  355,  356 

CSUPER,  79,  4,  345,  346,  355,  356,  445,  446 
,  455,  456 

CSUPER,  80,  49,  445,  446,  455,  456,  545,  546 
,  555,  556 

$  ****:M:************************************** 

$  GLOBAL  DUMMY  ELEMENT 

CHEXA,  999,  31,  111,  511,  551,  151,  116,  516 

,  556,  156 

PSOLID,  31,  41 

MAT1,  41,  3.0+6,  ,  0.0,  1.0-8 

$  GLOBAL  RESTRAINTS  AND  FORCES 


SPC, 

20, 

111, 

123456 

SPC, 

20, 

121, 

3 

SPC, 

20, 

131, 

3 

SPC, 

20, 

141, 

3 

SPC, 

20, 

151, 

3 

SPC, 

20, 

211, 

3 

SPC, 

20, 

221, 

3 

SPC, 

20, 

231, 

3 

SPC, 

20, 

241, 

3 

SPC, 

20, 

251, 

3 

SPC, 

20, 

311, 

3 

SPC, 

20, 

321, 

3 

SPC, 

20, 

331, 

3 

SPC, 

20, 

341, 

3 

SPC, 

20, 

351, 

3 

SPC, 

20, 

411, 

3 

SPC, 

20, 

421, 

3 

SPC, 

20, 

431, 

3 

SPC, 

20, 

441, 

3 

SPC, 

20, 

451, 

3 

SPC, 

20, 

551, 

3 

SPC, 

20, 

541, 

3 

SPC, 

20, 

531, 

3 

SPC, 

20, 

521, 

3 

SPC, 

20, 

511, 

23456 

PARAM,  POST,  -2 
PARAM,  AUTOS  PC,  YES 
FORCE,  10,  116,  ,  2.0, 

FORCE,  10,  126,  ,  4.0, 

FORCE,  10,  136,  ,  4.-:, 

FORCE,  10,  146,  ,  4.0, 

FORCE,  10,  156,  ,  2.0, 

FORCE,  10,  216,  ,  4.0, 

FORCE,  10,  226,  ,  8.0, 


1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 


FORCE,  10,  236,  ,  8.0,  ,  ,  1.0 


96 


FORCE, 

10, 

246, 

1 

8.0, 

y 

1.0 

FORCE, 

10, 

256, 

» 

4.0, 

y 

1.0 

FORCE, 

10, 

316, 

f 

4.0, 

y 

1.0 

FORCE, 

10, 

326, 

1 

8.0, 

y 

1.0 

FORCE, 

10, 

336, 

> 

8.0, 

y 

1.0 

FORCE, 

10, 

346, 

» 

8.0, 

y 

1.0 

FORCE, 

10, 

356, 

t 

4.0, 

y 

1.0 

FORCE, 

10, 

416, 

9 

4.0, 

9 

1.0 

FORCE, 

10, 

426, 

} 

8.0, 

y 

1.0 

FORCE, 

10, 

J 

8.0, 

y 

1.0 

FORCE, 

10, 

4**0, 

) 

8.0, 

y 

1.0 

FORCE, 

10, 

456, 

I 

4.0, 

y 

1.0 

FORCE, 

10, 

516, 

1 

2.0, 

y 

1.0 

FORCE, 

10, 

526, 

y 

4.0, 

y 

1.0 

FORCE, 

10, 

536, 

y 

4.0, 

y 

1.0 

FORCE, 

10, 

546, 

y 

4.0, 

y 

1.0 

FORCE, 

10, 

556, 

y 

2.0, 

y 

1.0 

ENDDATA 

97 


Appendix  B:  Computer  Model  Test  Runs.  Input  Data  and  Results 


Tests:  PORE1- 

Unit  length  =  1  micrometer  (10“®  m) 

80  cells,  4x4x5,  each  with  dimension  20x20x20  units 
Eo  =  29.3  x  106  psi&,  2.02  x  10u  Pascals 
Applied  force,  upward:  128  x  10“u  Newtons 
Evenly  distributed  pores  (one  per  cell) 

Simulated  material:  Porous  Cobolt  Oxide  (CoO) 


Data: 

Test  Name 

Avg  Disp 
li  IQiHl 

Calculated 
£  UOHEai 

%  Porosity 

Pore  Size 

Pore 1-1 

9.9 

2.0200 

0 

N/A 

Porel-2 

10.7 

1.8692 

6.4 

8  units 

Porel-3 

11.1 

1.8018 

9.1 

9  units 

Pore 1-4 

12.36 

1.6181 

16.6 

11  units 

Porel-5 

13.3 

1.5138 

21.6 

12  units 

Tests:  PORE1-  B 


Unit  length  =  1  micrometer  (10-6  m) 

80  cells,  4x4x5,  each  with  dimension  20x20x20  units 
Eo  =  33.2  x  106  psia,  2.29  x  1011  Pascals 
Applied  force,  upward:  128  x  10-u  Newtons 
Evenly  distributed  pores  (one  per  cell) 

Simulated  material:  Porous  Cobolt  Oxide  (CoO) 


Data: 

• 

Avg  Disp 

Test  Name 

ix  10-12) 

Porel-IB 

8.7336 

Porel-2B 

9.4574 

Porel-3B 

9.8110 

• 

Porel-4B 

10.9318 

Porel-5B 

11.7966 

Calculated 

E  (10u  Pa) 

%  Porosi 

fcZ  Pore  Size 

2.2900 

0 

N/A 

2.1147 

6.4 

8  units 

2.0385 

9.1 

9  units 

1.8295 

16.6 

11  units 

1.6954 

21.6 

12  units 
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Test 


name:  PORE2-1 

Unit  length  =  1  micrometer  (10**  m) 

#  80  cells,  4x4x5,  each  with  dimension  20x20x20  units 
Eo  =  29.3  x  10*  psia,  2.02  x  10u  Pascals 

8.88  %  porosity 

Pore  size:  cubes,  8  units  on  each  side 
Simulated  material:  Porous  Cobolt  Oxide  (CoO) 

Pores  follow  a  poison  distribution: 

^  Pores  per  cell:  012345678 

Number  of  cells:  23  22  21  10  3  1 

Data: 

listed  below  are  the  displacements  of  the  25  node  points  on  the  top  of 
the  model  (each  x  10  _u  units). 

#  1  X  ->  Y 


1.275 

1.243 

1.204 

1.205 

1.224 

1.221 

1.214 

1.196 

1.200 

1.222 

1.174 

1.200 

1.211 

1.203 

1.203 

1.154 

1.188 

1.206 

1.202 

1.203 

1.131 

1.160 

1.173 

1.191 

1.216 

Results: 

average  displacement  =  12.181  x  10-12  units 
calculated  value  for  E  =  1.642  x  10u  Pascals 

=  23.813  x  10*  Psia 


•  Test  name:  PORE2-2 

Unit  length  =  1  micrometer  (10“®  m) 

80  cells,  4x4x5,  each  with  dimension  20x20x20  units 
Eo  =  29.3  x  10*  psia,  2.02  x  10u  Pascals 
8.00  %  porosity 

Pore  size:  cubes,  8  units  on  each  side 

•  Simulated  material:  Porous  Cobolt  Oxide  (CoO) 

Pores  follow  a  poison  distribution: 

Pores  per  cell:  012345678 
Number  of  cells:  26  23  20  8  2  1 

Data: 

•  listed  below  are  the  displacements  of  the  25  node  points  on  the  top  of 
the  model  (each  x  10  _u  units). 

IX  ->  Y 


1.285 

1.237 

1.183 

1.173 

1.181 

1.221 

1.198 

1.167 

1.159 

1.170 

1.166 

1.177 

1.174 

1.154 

1.143 

1.143 

1.163 

1.166 

1.149 

1.137 

1.118 

1.132 

1.131 

1.134 

1.146 

Results: 

average  displacement  =  11.684  x  10" 12  units 
calculated  value  for  E  =  1.711  x  10u  Pascals 

=  24.8  x  10*  Psia 
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Test 


name:  PORE2-3 

Unit  length  =  1  micrometer  ( 10”6  m) 

•  80  cells,  4x4x5,  each  with  dimension  20x20x20  units 
Eo  =  29.3  x  10®  psia,  2.02  x  10u  Pascals 

10.08  %  porosity 

Pore  size:  cubes,  8  units  on  each  side 
Simulated  material:  Porous  Cobolt  Oxide  (CoO) 

Pores  follow  a  poison  distribution: 

®  Pores  per  cell:  012345678 

Number  of  cells:  19  21  22  12  5  1 

Data: 

listed  below  are  the  displacements  of  the  25  node  points  on  the  top  of 
the  model  (each  x  10  "u  units). 

•  IX  ->Y 


1.578 

1.393 

1.233 

1.203 

1.179 

1.392 

1.293 

1.234 

1.203 

1.175 

1.244 

1.250 

1.244 

1.207 

1.172 

1.294 

1.263 

1.217 

1.191 

1.167 

1.462 

1.254 

1.171 

1.159 

1.151 

Results: 

average  displacement  =  12.53  x  10"12  units 
calculated  value  for  E  =  1.596  x  10u  Pascals 

=  23.15  x  10®  Psia 


•  Test  name:  PORE2-4 

Unit  length  =  1  micrometer  (10-®  m) 

80  cells,  4x4x5,  each  with  dimension  20x20x20  units 
Eo  =  29.3  x  10®  psia,  2.02  x  10u  Pascals 
11.52  %  porosity 

Pore  size:  cubes,  8  units  on  each  side 

•  Simulated  material:  Porous  Cobolt  Oxide  (CoO) 

Pores  follow  a  poison  distribution: 

Pores  per  cell:  012345678 
Number  of  cells:  15  19  23  15  6  2 

Patau 

•  listed  below  are  the  displacements  of  the  25  node  points  on  the  top  of 
the  model  (each  x  10  *u  units). 

i  X  ->  Y 


1.690 

1.526 

1.454 

1.459 

1.335 

1.459 

1.380 

1.367 

1.334 

1.274 

1.274 

1.297 

1.300 

1.265 

1.229 

1.294 

1.275 

1.240 

1.217 

1.195 

1.434 

1.236 

1.161 

1.155 

1.149 

Results: 

average  displacement  =  1.320  x  10-12  units 
calculated  value  for  E  =  1.515  x  1011  Pascals 

=  21.975  x  10®  Psia 


101 


Test 


name:  PORE2-5 

Unit  length  =  1  micrometer  (l(Hm) 

80  cells,  4x4x5,  each  with  dimension  20x20x20  units 
Eo  =  29.3  x  106  psia,  2.02  x  1011  Pascals 
14.16  %  porosity 

Pore  size:  cubes,  8  units  on  each  side 
Simulated  material:  Porous  Cobolt  Oxide  (CoO) 

Pores  follow  a  poison  distribution: 

w  Pores  per  cell:  012345678 

Number  of  cells:  10  15  23  18  9  4  1 

Data: 

listed  below  are  the  displacements  of  the  25  node  points  on  the  top  of 
the  model  (each  x  10  _u  units). 

•  l  X  ->  Y 


1.713 

1.560 

1.514 

1.592 

1.532 

1.503 

1.434 

1.440 

1.446 

1.426 

1.424 

1.412 

1.381 

1.358 

1.338 

1.600 

1.437 

1.339 

1.316 

1.289 

1.660 

1.391 

1.284 

1.265 

1.241 

•  Results: 

average  displacement  =  14.933  x  10"12  units 
calculated  value  for  E  *  1.339  x  10u  Pascals 

=  19.42  x  10*  Psia 


0  Test  name:  PORE2-6 

Unit  length  a  1  micrometer  (10~®  m) 

80  cells,  4x4x5,  each  with  dimension  20x20x20  units 
Eo  =  29.3  x  10*  psia,  2.02  x  10u  Pascals 
17.12  X  porosity 

Pore  size:  cubes,  8  units  on  each  side 
9  Simulated  material:  Porous  Cobolt  Oxide  (Co O) 

Pores  follow  a  poison  distribution: 

Pores  per  cell:  012345678 
Number  of  cells:  6  11  21  20  13  6  2  1 

Patft? 

•  listed  below  are  the  displacements  of  the  25  node  points  on  the  top  of 

the  model  (each  x  10  _u  units), 
i  X  ->  Y 


1.712 

1.618 

1.633 

1.741 

1.660 

1.528 

1.526 

1.593 

1.613 

1.580 

1.462 

1.517 

1.546 

1.550 

1.538 

1.643 

1.531 

1.487 

1.503 

1.504 

1.707 

1.476 

1.408 

1.423 

1.427 

Results: 

average  displacement  =  15.571  x  10*12  units 
calculated  value  for  E  =  1.28  x  10u  Pascals 

=  18.63  x  10*  Psia 
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Test  name:  PORE2-7 

Unit  length  =  1  micrometer  (10“®  m) 

80  cells,  4x4x5,  each  with  dimension  20x20x20  units 
Eo  =  29.3  x  10*  psia,  2.02  x  1011  Pascals 
18.08  X  porosity 

Pore  size:  cubes,  8  units  on  each  side 
Simulated  material:  Porous  Cobolt  Oxide  (CoO) 

Pores  follow  a  poison  distribution: 

Pores  per  cell:  012345678 
Number  of  cells:  5  10  20  20  14  7  3  1 

Data: 

listed  below  are  the  displacements  of  the  25  node  points  on  the  top  of 
the  model  (each  x  10  *u  units), 
t  X  ->  Y 


1.730 

1.653 

1.686 

1.812 

1.746 

1.536 

1.547 

1.630 

1.671 

1.653 

1.456 

1.555 

1.599 

1.588 

1.598 

1.629 

1.584 

1.548 

1.559 

1.540 

1.693 

1.471 

1.488 

1.497 

1.448 

Results: 

average  displacement  =  15.79  x  10-12  units 
calculated  value  for  E  =  1.266  x  10u  Pascals 

=  18.36  x  10®  Psia 


Test  name:  PORE2-8 

Unit  length  a  1  micrometer  (10“®  m) 

80  cells,  4x4x5,  each  with  dimension  20x20x20  units 
Eo  =  29.3  x  10®  psia,  2.02  x  10u  Pascals 
19.04  X  porosity 

Pore  size:  cubes,  8  units  on  each  side 
Simulated  material:  Porous  Cobolt  Oxide  (CoO) 

Pores  follow  a  poison  distribution: 


Pores  per  cell:  0 
Number  of  cells:  4 


2  3  4 
19  20  15 


Data: 

listed  below  are  the  displacements  of  the  25  node  points  on  the  top  of 
the  model  (each  x  10  "u  units), 
i  X  ->  Y 


1.725 

1.796 

1.854 

1.832 

1.771 

1.534 

1.597 

1.684 

1.694 

1.681 

1.461 

1.561 

1.613 

1.612 

1.627 

1.633 

1.597 

1.568 

1.620 

1.664 

1.693 

1.483 

1.508 

1.660 

1.682 

Results: 

average  displacement  =  16.46 
calculated  value  for  E  =  1.215 

=  17.62 


10* 12  units 
10u  Pascals 
10®  Psia 
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Test  name:  PORE2-9 

Unit  length  =  1  micrometer  (10-6  m) 

•  80  cells,  4x4x5,  each  with  dimension  20x20x20  units 

Eo  =  29.3  x  10®  psia,  2.02  x  10u  Pascals 
20.56  X  porosity 

Pore  size:  cubes,  8  units  on  each  side 
Simulated  material:  Porous  Cobolt  Oxide  (CoO) 

Pores  follow  a  poison  distribution: 

^  Pores  per  cell:  012345678 

Number  of  cells:  3  8  17  20  16  9  4  2  1 


Data: 

listed  below  are  the  displacements  of  the  25  node 
model  (each  x  10  “u  units). 

X  ->  Y 


the 

i 


1.876 

1.911 

1.937 

2.001 

2.070 

1.730 

1.741 

1.777 

1.797 

1.790 

1.627 

1.686 

1.693 

1.668 

1.674 

1.709 

1.656 

1.609 

1.650 

1.685 

1.711 

1.490 

1.508 

1.655 

1.669 

Results: 

average  displacement  =  17.33 
calculated  value  for  E  =  1.154 

=  16.73 


x  10"12  units 
x  10u  Pascals 
x  10®  Psia 


points  on  the  top  of 


104 


Test 


name:  PORE2-2;2 

Unit  length  =  1  micrometer  (10“®  m) 

80  cells,  4x4x5,  each  with  dimension  20x20x20  units 
Eo  =  26.1  x  106  psia,  1.8  x  1011  Pascals 
8.00  X  porosity 

Pore  size:  cubes,  8  units  on  each  side 
Simulated  material:  Porous  Cobolt  Oxide  (CoO) 

Pores  follow  a  poison  distribution: 

Pores  per  cell:  012345678 
Number  of  cells:  26  23  20  8  2  1 


Data: 

listed  below  are  the  displacements  of  the  25  node  points  on  the  top  of 


the  model 

(each  x 

10  *u  units). 

• 

i  X  -> 

Y 

1.442 

1.388 

1.328 

1.316 

1.325 

1.370 

1.344 

1.309 

1.301 

1.312 

1.309 

1.322 

1.317 

1.295 

1.283 

1.283 

1.305 

1.309 

1.289 

1.276 

• 

1.255 

Results: 

1.270 

1.269 

1.273 

1.286 

average  displacement  =  13.112  x  10-12  units 
calculated  value  for  E  =  1.525  x  1011  Pascals 

=  22.12  x  10®  Psia 


Test  name:  PORE2-3;2 

Unit  length  =  1  micrometer  (10~®  m) 

80  cells,  4x4x5,  each  with  dimension  20x20x20  units 
Eo  =  26.1  x  10®  psia,  1.8  x  10u  Pascals 
10.08  X  porosity 

Pore  size:  cubes,  8  units  on  each  side 
Simulated  material:  Porous  Cobolt  Oxide  (CoO) 

Pores  follow  a  poison  distribution: 

Pores  per  cell:  012345678 
Number  of  cells:  19  21  22  12  5  1 

Data: 

listed  below  are  the  displacements  of  the  25  node  points  on  the  top  of 
the  model  (each  x  10  _u  units), 
i  X  ->  Y 


1.771 

1.563 

1.383 

1.350 

1.324 

1.562 

1.451 

1.385 

1.350 

1.319 

1.396 

1.402 

1.396 

1.355 

1.315 

1.452 

1.417 

1.366 

1.336 

1.309 

1.641 

1.407 

1.314 

1.301 

1.292 

Results: 

average  displacement  =  14.064  x  10"12  units 
calculated  value  for  E  =  1.422  x  10u  Pascals 

=  20.62  x  10®  Psia 
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Test  name:  PORE2-5;2 

Unit  length  =  1  micrometer  (10~*  m) 

#  80  cells,  4x4x5,  each  with  dimension  20x20x20  units 

Eo  =  26.1  x  10s  psia,  1.8  x  10u  Pascals 
14.16  %  porosity 

Pore  size:  cubes,  8  units  on  each  side 
Simulated  material:  Porous  Cobolt  Oxide  (CoO) 

Pores  follow  a  poison  distribution: 

®  Pores  per  cell:  012345678 

Number  of  cells:  10  15  23  28  9  4  1 


Data: 

listed  below  are  the  displacements  of  the  25  node  points  on  the  top  of 
the  model  (each  x  10  *u  units). 

•  1  X  ->  Y 


1.923 

1.751 

1.699 

1.787 

1.719 

1.687 

1.609 

1.616 

1.623 

1.600 

1.598 

1.584 

1.550 

1.525 

1.501 

1.796 

1.612 

1.503 

1.477 

1.447 

1.863 

1.561 

1.441 

1.419 

1.393 

Results: 

average  displacement  =  16.113  x  10"12  units 
calculated  value  for  E  =  1.241  x  10u  Pascals 

=  18.00  x  10*  Psia 


•  Teat  name:  PORE2-6;2 

Unit  length  *  1  micrometer  (10~*  m) 

80  cells,  4x4x5,  each  with  dimension  20x20x20  units 
Eo  =  26.1  x  10*  psia,  1.8  x  10u  Pascals 
17.12  X  porosity 

Pore  size:  cubes,  8  units  on  each  side 

•  Simulated  material:  Porous  Cobolt  Oxide  (CoO) 

Pores  follow  a  poison  distribution: 

Pores  per  cell:  012345678 
Number  of  cells:  6  11  21  20  13  6  2  1 

Data: 

•  listed  below  are  the  displacements  of  the  25  node  points  on  the  top  of 
the  model  (each  x  10  *u  units). 

i  X  ->  Y 


1.921 

1.815 

1.833 

1.954 

1.863 

1.714 

1.713 

1.787 

1.810 

1.773 

1.640 

1.703 

1.735 

1.739 

1.727 

1.844 

1.718 

1.669 

1.687 

1.687 

1.915 

1.656 

1.580 

1.579 

1.601 

Results: 

average  displacement  =  17.474  x  10-12  units 
calculated  value  for  E  =  1.145  x  10u  Pascals 

=  16.60  x  10«  Psia 
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Test 


name:  PORE2-8;2 

Unit  length  =  1  micrometer  (10*6  m) 

•  80  cells,  4x4x5,  each  with  dimension  20x20x20  units 

Eo  =  26.1  x  106  psia,  1.8  x  10u  Pascals 
19.04  %  porosity 

Pore  size:  cubes,  8  units  on  each  side 
Simulated  material:  Porous  Cobolt  Oxide  (CoO) 

Pores  follow  a  poison  distribution: 

^  Pores  per  cell:  012345678 

Number  of  cells:  4  9  19  20  15  8  4  1 


Data: 

listed  below  are  the  displacements  of  the  25  node 
the  model  (each  x  10  -u  units). 

•  l  X  ->  Y 


1.936 

2.016 

2.081 

1.056 

1.988 

1.721 

1.793 

1.889 

1.901 

1.886 

1.639 

1.752 

1.811 

1.809 

1.825 

1.832 

1.792 

1.760 

1.818 

1.868 

1.900 

1.664 

1.692 

1.863 

1.888 

Results: 

average  displacement  =  18.472  x  10-12  units 
calculated  value  for  E  =  1.0827  x  10u  Pascals 


=  15.70  x  10*  Psia 


points  on  the  top  of 
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Test  name:  PORE3-1 

Unit  length  =  .25  micrometer  (10-6  m) 

#  80  cells,  4x4x5,  each  with  dimension  20x20x20  units 
Eo  =  63.8  x  10®  psia,  4.4  x  10u  Pascals 

8.0  X  porosity 

Pore  size:  cubes,  8  units  on  each  side,  i.e.  2  micrometers 
Simulated  material:  Porous  Silicon  Carbide  (SiC) 

Pores  follow  a  poison  distribution: 

®  Pores  per  cell:  012345678 

Number  of  cells:  26  23  20  8  2  1 

Data: 

listed  below  are  the  displacements  of  the  25  node  points  on  the  top  of 
the  model  (each  x  10  *12  units). 

•  i  X  ->  Y 


5.918 

5.692 

5.452 

5.409 

5.448 

5.618 

5.505 

5.363 

5.335 

5.390 

5.376 

5.414 

5.395 

5.307 

5.265 

5.273 

5.348 

5.364 

5.285 

5.233 

5.163 

5.216 

5.206 

5.222 

5.278 

•  Results: 

average  displacement  =  5.379  x  10'12  units 
calculated  value  for  E  =  3.7181  x  10u  Pascals 

=  53.925  x  106  Psia 


Teat  name:  PORE3-2 

Unit  length  =  .25  micrometer  (10**  m) 

80  cells,  4x4x5,  each  with  dimension  20x20x20  units 
Eo  =  63.8  x  106  psia,  4.4  x  10u  Pascals 
10.08  %  porosity 

Pore  size:  cubes,  8  units  on  each  side,  i.e.  2  micrometers 
Simulated  material:  Porous  Silicon  Carbide  (SiC) 

Pores  follow  a  poison  distribution: 

Pores  per  cell:  012345678 
Number  of  cells:  19  21  22  12  5  1 

Data: 

listed  below  are  the  displacements  of  the  25  node  points  on  the  top  of 
the  model  (each  x  10  ~u  units). 

IX  ->  Y 


7.248 

6.410 

5.682 

5.551 

5.444 

6.405 

5.946 

5.667 

5.536 

5.423 

5.740 

5.743 

5,712 

5.553 

5.402 

5.974 

5.814 

5.598 

5.479 

5.372 

6.747 

5.792 

5,394 

5.335 

5.296 

Results: 

average  displacement  =  5.771  x  10-12  units 
calculated  value  for  E  =  3.466  x  1011  Pascals 

=  50.264  x  10®  Psia 
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Test  name:  PORE3-3 

Unit  length  =  .25  micrometer  (10-6  m) 

#  80  cells,  4x4x5,  each  with  dimension  20x20x20  units 

Eo  =  63.8  x  106  psia,  4.4  x  10u  Pascals 
14.16  %  porosity 

Pore  size:  cubes,  8  units  on  each  side,  i.e.  2  micrometers 
Simulated  material:  Porous  Silicon  Carbide  (SiC) 


• 

Pores 

Pores  per 

follow 

cell: 

a 

0 

poison 

1 

distribution: 

2  3  4 

5 

6 

7 

8 

Number  of 

Data: 

cells: 

10 

15 

23  28  9 

4 

1 

listed  below  are  the  displacements  of  the  25  node  points  on  the  top  of 
the  model  (each  x  10  ~12  units). 

•  i  X  ->  Y 


7.901 

7.197 

6.984 

7.353 

7.088 

6.945 

6.605 

6.627 

6.661 

6.591 

6.585 

6.499 

6.343 

6.252 

6.180 

7.383 

6.618 

6.165 

6.061 

5.952 

7.660 

6.429 

5.930 

5.839 

5.734 

•  Results: 

average  displacement  =  6.623  x  10-12  units 
calculated  value  for  E  =  3.019  x  1011  Pascals 

=  43.8  x  10«  Psia 


£  Test  name:  PORE3-4 

Unit  length  =  .25  micrometer  (10“*  m) 

80  cells,  4x4x5,  each  with  dimension  20x20x20  units 
Eo  =  63.8  x  106  psia,  4.4  x  10u  Pascals 
20.56  %  porosity 

Pore  size:  cubes,  8  units  on  each  side,  i.e.  2  micrometers 
0  Simulated  material:  Porous  Silicon  Carbide  (SiC) 

Pores  follow  a  poison  distribution: 

Pores  per  cell:  012345678 
Number  of  cells:  3  8  17  20  16  9  4  2  1 

Data: 

0  listed  below  are  the  displacements  of  the  25  node  points  on  the  top  of 

the  model  (each  x  10  _12  units). 

IX  ->  Y 


8.691 

8.856 

8.973 

9.266 

9.594 

7.997 

8.037 

8.193 

8.294 

8.291 

7.525 

7.762 

7.783 

7.681 

7.732 

7.884 

7.638 

7.423 

7.614 

7.784 

7.884 

6.890 

6.974 

7.643 

7.705 

Results: 

average  displacement  =  8.005  x  10-12  units 
calculated  value  for  E  =  2.498  x  10u  Pascals 

=  36.24  x  106  Psia 
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"est  name:  PORB3-5 

Unit  length  =  .25  micrometer  (10*6  m) 
f  80  cells,  4x4x5,  each  with  dimension  20x20x20  units 

Eo  =  63.8  x  106  psia,  4.4  x  1011  Pascals 
0.7  X  porosity 

Pore  size:  cubes,  8  units  on  each  side,  i.e.  2  micrometers 
Simulated  material:  Porous  Silicon  Carbide  (SiC) 

Pores  follow  a  poison  distribution: 

®  Pores  per  cell:  012345678 

Number  of  cells:  72  7  1 

Data: 

listed  below  are  the  displacements  of  the  25  node  points  on  the  top  of 
the  model  (each  x  10  ~12  units). 

•  i  X  ->  Y 


4.572 

4.562 

4.556 

4.563 

4.568 

4.590 

4.577 

4.567 

4.565 

4.561 

4.603 

4.617 

4.602 

4.566 

4.555 

4.621 

4.633 

4.616 

4.577 

4.563 

4.639 

4.620 

4.602 

4.591 

4.576 

•  Results: 

average  displacement  =  4.586  x  10-12  units 
calculated  value  for  E  =  4.361  x  1011  Pascals 

=  63.244  x  10«  P~: 


Test  name:  PORE3-6 

Unit  length  =  .25  micrometer  (10-6  m) 

80  cells,  4x4x5,  each  with  dimension  20x20x20  units 
Eo  =  63.8  x  106  psia,  4.4  x  10u  Pascals 
0.96  X  porosity 

Pore  size:  cubes,  8  units  on  each  side,  i.e.  2  micrometers 
Simulated  material:  Porous  Silicon  Carbide  (SiC) 

Pores  follow  a  poison  distribution: 

Pores  per  cell:  012345678 
Number  of  cells:  69  10  1 


Data: 

listed  below  are  the  displacements  of  the  25  node  points  on  the  top  of 


the  model 
t  X  -> 

(each  x 
Y 

10  *12  units). 

4.558 

4.558 

4.562 

4.575 

4.587 

4.581 

4.577 

4.577 

4.582 

4.584 

4.599 

4.621 

4.616 

4.588 

4.583 

4.622 

4.643 

4.635 

4.605 

4.600 

4.645 

Results: 

4.635 

4.626 

4,625 

4.618 

average  displacement  =  4.600  x  10-12  units 
calculated  value  for  E  =  4.3477  x  1011  Pascals 

=  63.056  x  106  Psia 
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Test 


name:  PORE3-7 

Unit  length  =  .25  micrometer  (10*®  m) 

#  80  cells,  4x4x5,  each  with  dimension  20x20x20  units 
Eo  =  63.8  x  106  psia,  4.4  x  10u  Pascals 

1.2  %  porosity 

Pore  size:  cubes,  8  units  on  each  side,  i.e.  2  micrometers 
Simulated  material:  Porous  Silicon  Carbide  (SiC) 

Pores  follow  a  poison  distribution: 

®  Pores  per  cell:  012345678 

Number  of  cells:  67  11  2 

Data: 

listed  below  are  the  displacements  of  the  25  node  points  on  the  top  of 
the  model  (each  x  10  -12  units). 

•  i  X  ->  Y 


4.576 

4.573 

4.574 

4.586 

4.597 

4.594 

4.589 

4.588 

4.591 

4.592 

4.648 

4.647 

4.621 

4.592 

4.588 

4.669 

4.664 

4.652 

4.621 

4.592 

4.643 

4.630 

4.662 

4.658 

4.604 

•  Results: 

average  displacement  =  4.614  x  10"12  units 
calculated  value  for  E  =  4.335  x  1011  Pascals 

=  62.866  x  106  Psia 


0  Test  name:  PORE 3-8 

Unit  length  =  .25  micrometer  (10-6  m) 

80  cells,  4x4x5,  each  with  dimension  20x20x20  units 
Eo  =  63.8  x  106  psia,  4.4  x  1011  Pascals 
2.16  %  porosity 

Pore  size:  cubes,  8  units  on  each  side,  i.e.  2  micrometers 
0  Simulated  material:  Porous  Silicon  Carbide  (SiC) 

Pores  follow  a  poison  distribution: 

Pores  per  cell:  012345678 
Number  of  cells:  59  16  4  1 

Data: 

•  listed  below  are  the  displacements  of  the  25  node  points  on  the  top  of 


the  model 
i  X  -> 

(each  x 
Y 

10  _12  units). 

4.941 

4.822 

4.705 

4.664 

4.645 

4.835 

4.755 

4.704 

4.672 

4.625 

4.769 

4.738 

4.707 

4.662 

4.618 

4.732 

4.718 

4.699 

4.663 

4.627 

4.670 

Results: 

4.654 

4.692 

4.691 

4.629 

^  m  n 

average  displacement  =  4.7055  x  10-12  units 
calculated  value  for  E  =  4.250  x  10u  Pascals 

=  61.64  x  106  Psia 
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Test 


none:  PORE3-9 

Unit  length  =  .25  micrometer  (1(H  m) 

•  80  cells,  4x4x5,  each  with  dimension  20x20x20  units 

Eo  =  63.8  x  106  psia,  4.4  x  10u  Pascals 
2.688  X  porosity 

Pore  size:  cubes,  8  units  on  each  side,  i.e.  2  micrometers 
Simulated  material:  Porous  Silicon  Carbide  (SiC) 

Pores  follow  a  poison  distribution: 

^  Pores  per  cell:  012345678 

Number  of  cells:  50  10  8  2 


Data: 

listed  below  are  the  displacements  of  the  25  node  points  on  the  top  of 
the  model  (each  x  10  _12  units). 

•  1  X  ->  Y 


4.925 

4.835 

4.742 

4.721 

4.725 

4.843 

4.799 

4.780 

4.778 

4.772 

4.802 

4.794 

4.793 

4.842 

4.909 

4.762 

4.805 

4.813 

4.835 

4.910 

4.683 

4.783 

4.834 

4.771 

4.758 

Results: 

average  displacement  =  4.8006  x  10-12  units 


calculated  value  for  E  =  4.166  x  10u  Pascals 

=  60.422  x  106  Psia 
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Test  name:  PORE3-l;2 

Unit  length  =  .25  micrometer  (10-6  m) 

80  cells,  4x4x5,  each  with  dimension  20x20x20  units 
Eo  =  61.64  x  106  psia,  4.25  x  10u  Pascals 
8.0  %  porosity 

Pore  size:  cubes,  8  units  on  each  side,  i.e.  2  micrometers 
Simulated  material:  Porous  Silicon  Carbide  (SiC) 

Pores  follow  a  poison  distribution: 

Pores  per  cell:  012345678 
Number  of  cells:  26  23  20  8  2  1 


Data: 

listed  below  are  the  displacements  of  the  25  node  points  on  the  top  of 
the  model  (each  x  10  ~n  units). 

IX  ->  Y 


6.127 

5.893 

5.645 

5.600 

5.640 

5.817 

5.699 

5.553 

5.523 

5.580 

5.565 

5.605 

5.586 

5.494 

5.451 

5.459 

5.537 

5.553 

5.471 

5.418 

5.345 

5.400 

5.390 

5.406 

5.464 

Results: 

average  displacement  =  5.569  x  10"12  units 
calculated  value  for  E  =  3.590  x  1011  Pascals 

=  52.08  x  10®  Psia 


•  Test  name:  PORE3-2;2 

Unit  length  =  .25  micrometer  (10“®  m) 

80  cells,  4x4x5,  each  with  dimension  20x20x20  units 
Eo  =  61.64  x  106  psia,  4.25  x  10u  Pascals 
10.08  %  porosity 

Pore  size:  cubes,  8  units  on  each  side,  i.e.  2  micrometers 

•  Simulated  material:  Porous  Silicon  Carbide  (SiC) 

Pores  follow  a  poison  distribution: 

Pores  per  cell:  012345678 
Number  of  cells:  19  21  22  12  5  1 

Data: 

•  listed  below  are  the  displacements  of  the  25  node  points  on  the  top  of 
the  model  (each  x  10  ~u  units). 

IX  ->  Y 


7.504 

6.636 

5.883 

5.747 

5.636 

6.631 

6.156 

5.870 

5.732 

5.615 

5.943 

5.946 

5.914 

5.749 

5.593 

6.184 

6,019 

5.796 

5.673 

5.562 

6.985 

5.997 

5.584 

5.524 

5.483 

Results: 

average  displacement  =  5.974  x  10-12  units 
calculated  value  for  E  =  3.3475  x  10u  Pascals 

=  48.55  x  10®  Psia 
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Test  name:  PORE3-3;2 

Unit  length  =  .25  micrometer  ( 10-6  m) 

#  80  cells,  4x4x5,  each  with  dimension  20x20x20  units 
Eo  =  61.64  x  106  psia,  4.25  x  1011  Pascals 

14.16  %  porosity 

Pore  size:  cubes,  8  units  on  each  side,  i.e.  2  micrometers 
Simulated  material:  Porous  Silicon  Carbide  (SiC) 

Pores  follow  a  poison  distribution: 

®  Pores  per  cell:  012345678 

Number  of  cells:  10  15  23  18  9  4  1 

Data: 

listed  below  are  the  displacements  of  the  25  node  points  on  the  top  of 
the  model  (each  x  10  ~u  units). 

•  i  X  ->  Y 


8.180 

7.451 

7.231 

7.612 

7.339 

7.190 

6.838 

6.861 

6.897 

6.823 

6.817 

6.729 

6.567 

6.473 

6.398 

7.643 

6.852 

6.382 

6.275 

6.162 

7.930 

6.655 

6.139 

6.045 

5.936 

Results: 

average  displacement  =  6.857  x  10'12  units 
calculated  value  for  E  =  2.917  x  10u  Pascals 

«  42.3  x  10«  Psia 


•  Teat  name:  PORE3-4;2 

Unit  length  =  .25  micrometer  (10“®  m) 

80  cells,  4x4x5,  each  with  dimension  20x20x20  units 
Eo  =  61.64  x  106  psia,  4.25  x  10u  Pascals 
20.56  X  porosity 

Pore  size:  cubes,  8  units  on  each  side,  i.e.  2  micrometers 

•  Simulated  material:  Porous  Silicon  Carbide  (SiC) 

Pores  follow  a  poison  distribution: 

Pores  per  cell:  012345678 
Number  of  cells:  3  8  17  20  16  9  4  2  1 

Data: 

•  listed  below  are  the  displacements  of  the  25  node  points  on  the  top  of 
the  model  (each  x  10  ~u  units). 

IX  ->  Y 


8.997 

9.168 

9.289 

9.593 

9.933 

8.279 

8.321 

8.482 

8.587 

8.583 

7.790 

8.036 

8.058 

7.953 

8.005 

8.163 

7.907 

7.685 

7.883 

8.059 

8.162 

7.133 

7.221 

7.913 

7.977 

Results: 

average  displacement  =  8.287  x  10-12  units 
calculated  value  for  E  =  2.413  x  1011  Pascals 

=  35.00  x  106  Psia 
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Test  name:  PORE3-11 

Unit  length  =  .25  micrometer  (10**  m) 

^  80  cells,  4x4x5,  each  with  dimension  20x20x20  units 

Eo  =  61.64  x  106  psia,  4.25  x  10u  Pascals 
0.7  %  porosity 

Pore  size:  cubes,  8  units  on  each  side,  i.e.  2  micrometers 
Simulated  material:  Porous  Silicon  Carbide  (SiC) 

Pores  follow  a  poison  distribution: 

®  Pores  per  cell:  012345678 

Number  of  cells:  77  7  1 

Data: 

listed  below  are  the  displacements  of  the  25  node  points  on  the  top  of 
the  model  (each  x  10  -12  units). 

•  1  X  ->  Y 


4.733 

4.723 

4.717 

4.723 

4.729 

4.752 

4.739 

4.728 

4.726 

4.722 

4.765 

4.780 

4.764 

4.728 

4.716 

4.784 

4.796 

4.778 

4.738 

4.724 

4.802 

4.783 

4.764 

4.753 

4.737 

•  Results: 

average  displacement  =  4.748  x  10”12  units 
calculated  value  for  E  =  4.212  x  10u  Pascals 

=  61.088  x  106  Psia 


Teat  name:  PORE 3- 12 

Unit  length  =  .25  micrometer  (10*6  m) 

80  cells,  4x4x5,  each  with  dimension  20x20x20  units 
Eo  =  61.64  x  106  psia,  4.25  x  10u  Pascals 
1.2  X  porosity 

Pore  size:  cubes,  8  units  on  each  side,  i.e.  2  micrometers 
Simulated  material:  Porous  Silicon  Carbide  (SiC) 

Pores  follow  a  poison  distribution: 

Pores  per  cell:  012345678 
Number  of  cells:  67  11  2 


Data: 

listed  below  are  the  displacements  of  the  25  node  pointB  on  the  top  of 
the  model  (each  x  10  _12  units), 
i  X  ->  Y 


4.738 

4.734 

4.736 

4.749 

4.759 

4.756 

4,751 

4.750 

4.752 

4.754 

4.812 

4.811 

4.784 

4.754 

4.749 

4.834 

4.828 

4.816 

4.784 

4.754 

4.807 

4.794 

4.827 

4.823 

4.766 

Results: 

average  displacement  =  4.777  x  10-12  units 
calculated  value  for  E  =  4.187  x  10u  Pascals 

=  60.72  x  10*  Psia 
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Test 


name:  PORB3-13 

Unit  length  =  .25  micrometer  (1(H  m) 

•  80  cells,  4x4x5,  each  with  dimension  20x20x20  units 

Eo  =  61.64  x  10®  psia,  4.25  x  1011  Pascals 
2.16  X  porosity 

Pore  size:  cubes,  8  units  on  each  side,  i.e.  2  micrometers 
Simulated  material:  Porous  Silicon  Carbide  (SiC) 

Pores  follow  a  poison  distribution: 

®  Pores  per  cell:  012345678 

Number  of  cells:  59  16  4  1 

Data: 

listed  below  are  the  displacements  of  the  25  node  points  on  the  top  of 


the  model 

(each  x 

10  _12  units). 

• 

1  X  -> 

Y 

5.115 

4.992 

4.871 

4.829 

4.809 

5.005 

4.923 

4.870 

4.837 

4.788 

4.937 

4.906 

4.873 

4.826 

4.781 

4.899 

4.884 

4.864 

4.828 

4.790 

• 

4.835 

Results: 

4.819 

4.858 

4.856 

4.792 

average  displacement  =  4.873  x  10-12  units 
calculated  value  for  E  =  4.104  x  1011  Pascals 

=  59.52  x  10®  Psia 
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Test  name:  P0RE4-1 

Unit  length  =  1  micrometer  (10*6  m) 

#  80  cells,  4x4x5,  each  with  dimension  20x20x20  units 
Eo  =  46.7  x  106  psia,  3.22  x  10u  Pascals 

8.00  %  porosity 

Pore  size:  cubes,  8  units  on  each  side 
Simulated  material:  Porous  Aluminium  Nitride  (AIN) 

Pores  follow  a  poison  distribution: 

®  Pores  per  cell:  012345678 

Number  of  cells:  26  23  20  8  2  1 

Data: 

listed  below  are  the  displacements  of  the  25  node  points  on  the  top  of 
the  model  (each  x  10  _12  units). 

•  i  X  ->  Y 


7.643 

7.538 

7.414 

7.509 

7.667 

7.264 

7.304 

7.297 

7.406 

7.596 

6.968 

7.201 

7.342 

7.363 

7.434 

6.883 

7.150 

7.319 

7.345 

7.405 

6.776 

7.002 

7.128 

7.278 

7.484 

Results: 

average  displacement  =  7.3086  x  10" 12  units 
calculated  value  for  E  =  2.7365  x  10n  Pascals 

=  39.688  x  10«  Psia 


0  Teat  name:  PORE4-2 

Unit  length  =  1  micrometer  (10-®  m) 

80  cells,  4x4x5,  each  with  dimension  20x20x20  units 
Eo  =  46.7  x  106  psia,  3.22  x  10u  Pascals 
10.08  X  porosity 

Pore  size:  cubes,  8  units  on  each  side 

•  Simulated  material:  Porous  Aluminium  Nitride  (AIN) 

Pores  follow  a  poison  distribution: 

Pores  per  cell:  012345678 
Number  of  cells:  19  21  22  12  5  1 

Data: 

•  listed  below  are  the  displacements  of  the  25  node  points  on  the  top  of 
the  model  (each  x  10  _12  units). 

IX  ->  Y 


8.444 

7.614 

6.904 

7.026 

7.317 

7.560 

7.263 

7.436 

7.946 

8.296 

6.930 

7.292 

7.898 

8.423 

8.738 

7.556 

7.715 

7.815 

8.166 

8.586 

8.936 

7.971 

7.817 

8.185 

8.615 

Results: 

average  displacement  =  7.858  x  10-12  units 
calculated  value  for  E  =  2.545  x  10u  Pascals 

=  36.91  x  10®  Psia 


Test 


name:  PORE4-3 

Unit  length  =  1  micrometer  (lO^m) 

80  cells,  4x4x5,  each  with  dimension  20x20x20  units 
Eo  =  46.7  x  106  psia,  3.22  x  1011  Pascals 
11.52  X  porosity 

Pore  size:  cubes,  8  units  on  each  side 
Simulated  material:  Porous  Aluminium  Nitride  (AIN) 

Pores  follow  a  poison  distribution: 

•  Pores  per  cell:  012345678 

Number  of  cells:  15  19  23  15  6  2 

Data: 

listed  below  are  the  displacements  of  the  25  node  points  on  the  top  of 
the  model  (each  x  10  ~n  units). 

•  i  X  ->  Y 


8.110 

8.051 

8.389 

8.672 

8.161 

7.676 

7.848 

8.185 

8.253 

8.136 

7.311 

7.752 

8.141 

8.229 

8.215 

7.704 

8.160 

8.350 

8.501 

9.020 

8.905 

8.825 

8.300 

8.855 

9.724 

#  Results: 

average  displacement  =  8.181  x  10" 12  units 
calculated  value  for  E  =  2.444  x  1011  Pascals 

=  35.455  x  106  Psia 


Test  name:  PORE4-4 

Unit  length  =  1  micrometer  (10“*  m) 

80  cells,  4x4x5,  each  with  dimension  20x20x20  units 
Eo  =  46.7  x  106  psia,  3.22  x  10u  Pascals 
14.16  X  porosity 

Pore  size:  cubes,  8  units  on  each  side 
Simulated  material:  Porous  Aluminium  Nitride  (AIN) 

Pores  follow  a  poison  distribution: 

Pores  per  cell:  012345678 
Number  of  cells:  10  15  23  18  9  4  1 

Data: 

listed  below  are  the  displacements  of  the  25  node  points  on  the  top  of 
the  model  (each  x  10  units). 

IX  ->  Y 


8.395 

8.452 

9.005 

9.721 

9.614 

8.251 

8.380 

8.739 

9.044 

9.223 

8.567 

8.619 

8.746 

8.928 

8.874 

9.700 

8.936 

8.797 

9.053 

9.372 

10.325 

8.935 

8.486 

9.192 

10.035 

Results: 

average  displacement  =  9.016  x  10"12  units 
calculated  value  for  E  =  2.218  x  1011  Pascals 

=  32.173  x  10«  Psia 
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Test  name:  PORE4-5 

Unit  length  =  1  micrometer  ( 10“*  m) 

%  80  cells,  4x4x5,  each  with  dimension  20x20x20  units 

Eo  =  46.7  x  106  psia,  3.22  x  10u  Pascals 
17.12  %  porosity 

Pore  size:  cubes,  8  units  on  each  side 
Simulated  material:  Porous  Aluminium  Nitride  (AIN) 

Pores  follow  a  poison  distribution: 

®  Pores  per  cell:  012345678 

Number  of  cells:  6  11  21  20  13  6  2  1 

Data: 

listed  below  are  the  displacements  of  the  25  node  points  on  the  top  of 
the  model  (each  x  10  “ 12  units). 


X  -> 

Y 

9.425 

9.462 

9.837 

10.060 

9.314 

8.987 

9.274 

10.046 

10.009 

9.340 

9.092 

9.410 

9.862 

9.909 

9.667 

10.269 

9.614 

9.384 

9.858 

10.512 

10.772 

9.383 

9.071 

9.974 

10.965 

Results: 

average  displacement  =  9.740  x  10-12  units 
calculated  value  for  E  =  2.054  x  1011  Pascals 

=  29.78  x  106  p8ia 


^  Teat  name:  PORE4-6 

Unit  length  =  1  micrometer  (10-6  m) 

80  cells,  4x4x5,  each  with  dimension  20x20x20  units 
Eo  =  46.7  x  10*  psia,  3.22  x  1011  Pascals 
19.04  X  porosity 

Pore  size:  cubes,  8  units  on  each  side 
0  Simulated  material:  Porous  Aluminium  Nitride  (AIN) 

Pores  follow  a  poison  distribution: 

Pores  per  ceil:  012345678 
Number  of  cells:  4  9  19  20  15  8  4  1 

Data: 

•  listed  below  are  the  displacements  of  the  25  node  points  on  the  top  of 

the  model  (each  x  10  ~u  units). 


X  -> 

Y 

8.944 

9.913 

11.361 

11.319 

11.084 

9.175 

9.837 

10.328 

10.559 

10.652 

9.876 

9.870 

10.271 

10.418 

10.482 

10.029 

9.904 

10.162 

10.832 

11.113 

10.520 

9.309 

9.610 

10.997 

11.674 

Results: 

average  displacement  =  10.330  x  10-12  units 
calculated  value  for  E  =  1.936  x  1011  Pascals 

=  28.08  x  106  Psia 
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Test  name:  PORE4-7 

Unit  length  =  1  micrometer  (10-*  m) 

80  cells,  4x4x5,  each  with  dimension  20x20x20  units 
Eo  =  46.7  x  10*  psia,  3.22  x  1011  Pascals 
20.56  X  porosity 

Pore  size:  cubes,  8  units  on  each  side 
Simulated  material:  Porous  Aluminium  Nitride  (AIN) 

Pores  follow  a  poison  distribution: 

Pores  per  cell:  012345678 
Number  of  cells:  3  8  17  20  16  9  4  2  1 

Data: 

listed  below  are  the  displacements  of  the  25  node  points  on  the  top 


the  model 

(each  x 

10  *u  units). 

• 

i  X  -> 

Y 

1.183 

1.206 

1.222 

1.262 

1.306 

1.089 

1.096 

1.118 

1.132 

1.129 

1.025 

1.060 

1.064 

1.049 

1.053 

1.075 

1.042 

1.013 

1.039 

1.061 

1.076 

0.938 

0.949 

1.042 

1.051 

• 

Results: 

x  10-12  units 

average  displacement  = 

10.013 

calculated  value  for  E 

=  1.8326 

x  1011  Pascals 

=  26.58 

x  106  Psia 
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Test  name:  PORE5-1 

Unit  length  =  1  micrometer  (10-®  m) 

80  cells,  4x4x5,  each  with  dimension  20x20x20  units 
Eo  =  46.7  x  106  psia,  3.22  x  10u  Pascals 
15.4  %  porosity 

Pore  size:  cubes,  10  units  on  each  side 
Simulated  material:  Porous  Aluminium  Nitride  (AIN) 

Pores  follow  a  poison  distribution: 

Pores  per  cell:  012345678 
Number  of  cells:  26  23  20  8  2  1 

Data: 

listed  below  are  the  displacements  of  the  25  node  points  on  the  top  of 
the  model  (each  x  10  ~12  units), 
i  X  ->  Y 


10.364 

9.881 

9.477 

9.787 

10.241 

8.980 

8.982 

8.956 

9.267 

9.802 

7.952 

8.476 

8.808 

8.912 

9.142 

7.591 

8.221 

8.617 

8.713 

8.887 

7.183 

7.737 

8.065 

8.435 

8.929 

Results: 

average  displacement  =  8.856  x  10"12  units 
calculated  value  for  E  =  2.258  x  10u  Pascals 

=  32.753  x  10®  Psia 


Test  name:  PORE5-2 

Unit  length  =  1  micrometer  (10“®  m) 

80  cells,  4x4x5,  each  with  dimension  20x20x20  units 
Eo  =  46.7  x  10®  psia,  3.22  x  10u  Pascals 
3.6  X  porosity 

Pore  size:  cubes,  10  units  on  each  side 
Simulated  material:  Porous  Aluminium  Nitride  (AIN) 

Pores  follow  a  poison  distribution: 

Pores  per  cell:  012345678 
Number  of  cells:  61  15  4 

Data: 

listed  below  are  the  displacements  of  the  25  node  points  on  the  top  of 
the  model  (each  x  10  *12  units), 
i  X  ->  Y 


6.310 

6.352 

6.398 

6.396 

6.389 

6.500 

6.449 

6.484 

6.582 

6.604 

6.569 

6.507 

6.512 

6.595 

6.619 

6.466 

6.546 

6.543 

6.480 

6.482 

6.445 

6.600 

6.616 

6.493 

6.483 

Results: 

average  displacement  =  6.497  x  10' 12  units 
calculated  value  for  E  =  3.078  x  1011  Pascals 

=  44.647  x  10®  Psia 
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Test  name:  PORE5-3 

Unit  length  =  1  micrometer  (10-®  m) 

•  80  cells,  4x4x5,  each  with  dimension  20x20x20  units 
Eo  =  46.7  x  106  psia,  3.22  x  10u  Pascals 

5.156  %  porosity 

Pore  size;  cubes,  10  units  on  each  side 
Simulated  material:  Porous  Aluminium  Nitride  (AIN) 

Pores  follow  a  poison  distribution: 

®  Pores  per  cell:  012345678 

Number  of  cells:  55  18  6  1 

Data: 

listed  below  are  the  displacements  of  the  25  node  points  on  the  top  of 
the  model  (each  x  10  _12  units). 

•  i  X  ->  Y 


7.643 

7.332 

6.985 

6.675 

6.443 

7.488 

7.151 

6.885 

6.752 

6.576 

7.156 

6.908 

6.756 

6.695 

6.564 

6.754 

6.723 

6.646 

6.503 

6.395 

6.515 

6.584 

6.563 

6.406 

6.321 

•  Results: 

average  displacement  =  6.777  x  10*12  units 
calculated  value  for  E  =  2.951  x  10*1  Pascals 

=  42.825  x  106  Psia 


£  Test  naae:  PORE5-4 

Unit  length  a  1  micrometer  (10"*  m) 

80  cells,  4x4x5,  each  with  dimension  20x20x20  units 
Eo  =  46.7  x  106  psia,  3.22  x  10u  Pascals 
6.56  %  porosity 

Pore  size:  cubes,  10  units  on  each  side 
9  Simulated  material:  Porous  Aluminium  Nitride  (AIN) 

Pores  follow  a  poison  distribution: 

Pores  per  cell:  012345678 
Number  of  cells:  50  20  8  2 

Data: 

•  listed  below  are  the  displacements  of  the  25  node  points  on  the  top  of 

the  model  (each  x  10  _12  units). 

IX  ->  Y 


7.644 

7.494 

7.270 

7.041 

6.895 

7.445 

7.259 

7.097 

7.059 

6.993 

7.088 

6.955 

7.026 

7.075 

6.910 

6.673 

6.751 

6.897 

6.856 

6.703 

6.419 

6.617 

6.677 

6.604 

6.636 

Results: 

average 

displacement  = 

6.963 

x  10'12  units 

calculated  value  for  E 

=  2.872 

x  10u  Pascals 

=  41.656  x  10®  Psia 
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Test  name:  PORE5-5 

Unit  length  =  1  micrometer  (10-6  m) 

80  cells,  4x4x5,  each  with  dimension  20x20x20  units 
Eo  =  46.7  x  106  psia,  3.22  x  10u  Pascals 
10.47  %  porosity 

Pore  size:  cubes,  10  units  on  each  side 
Simulated  material:  Porous  Aluminium  Nitride  (AIN) 

Pores  follow  a  poison  distribution: 

Pores  per  cell:  012345678 
Number  of  cells:  38  23  14  4  1 

Data: 

listed  below  are  the  displacements  of  the  25  node  points  on  the  top  of 
the  model  (each  x  10  units). 

IX  ->  Y 


9.622 

8.709 

8.042 

8.789 

10.420 

8.506 

8.061 

7.756 

8.045 

8.714 

7.336 

7.327 

7.455 

7.500 

7.449 

7.884 

7.303 

7.007 

7.108 

7.133 

7.721 

6.621 

6.833 

7.135 

Results: 

average  displacement  =  7.392  x  10”12  units 
calculated  value  for  E  =  2.706  x  10u  Pascals 

=  39.242  x  10*  Psia 


f  Teat  name:  PORE5-6 

Unit  length  =  1  micrometer  (10“®  m) 

80  cells,  4x4x5,  each  with  dimension  20x20x20  units 
Eo  =  46.7  x  106  psia,  3.22  x  1011  Pascals 
14.06  X  porosity 

Pore  size:  cubes,  10  units  on  each  side 

•  Simulated  material:  Porous  Aluminium  Nitride  (AIN) 

Pores  follow  a  poison  distribution: 

Pores  per  cell:  012345678 
Number  of  cells:  29  23  19  7  2 

Data: 

•  listed  below  are  the  displacements  of  the  25  node  points  on  the  top  of 
the  model  ( each  x  10  -*12  units). 

i  X  ->  Y 


10.477 

9.249 

8.664 

10,741 

9.105 

8.559 

8.281 

9.100 

10.300 

7.758 

7.845 

8.013 

8.034 

8.004 

8.405 

7.879 

7.626 

7.765 

7.837 

8.383 

7.316 

7.569 

7.915 

Results: 

average  displacement  =  8.277  x  10-12  units 
calculated  value  for  E  =  2.416  x  1011  Pascals 

=  35.046  x  10®  Psia 
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Test  name:  PORE5-7 

Unit  length  =  1  micrometer  (10*6  m) 

•  80  cells,  4x4x5,  each  with  dimension  20x20x20  units 
Eo  =  46.7  x  106  psia,  3.22  x  1011  Pascals 

19.06  %  porosity 

Pore  size:  cubes,  10  units  on  each  side 
Simulated  material:  Porous  Aluminium  Nitride  (AIN) 

Pores  follow  a  poison  distribution: 

®  Pores  per  cell:  012345678 

Number  of  cells:  20  21  22  12  4  1 

Data: 

listed  below  are  the  displacements  of  the  25  node  points  on  the  top  of 
the  model  (each  x  10  *u  units). 

•  i  X  ->  Y 


1.120 

1.049 

1.055 

1.321 

0.991 

0.972 

0.987 

1.111 

1.246 

0.870 

0.947 

0.999 

0.972 

1.004 

1.041 

1.032 

0.944 

0.931 

0.967 

1.059 

0.862 

0.889 

0.935 

Results: 

average  displacement  =  10.132  x  10~n  units 
calculated  value  for  E  =  1.974  x  10u  Pascals 

=  28.638  x  10s  Psia 


#  Test  name:  PORE5-8 

Unit  length  =  1  micrometer  (10-6  m) 

80  cells,  4x4x5,  each  with  dimension  20x20x20  units 
Eo  =  46.7  x  106  psia,  3.22  x  10u  Pascals 
16.875  %  porosity 

Pore  size:  cubes,  10  units  on  each  side 

#  Simulated  material:  Porous  Aluminium  Nitride  (AIN) 

Pores  follow  a  poison  distribution: 

Pores  per  cell:  012345678 
Number  of  cells:  24  22  21  9  3  1 

Data: 

#  listed  below  are  the  displacements  of  the  25  node  points  on  the  top  of 
the  model  (each  x  10  _1Z  units). 

IX  ->  Y 


11.901 

10.875 

10.560 

12.969 

10.013 

9.684 

9.749 

10.847 

12.039 

8.166 

8.437 

8.863 

9.092 

9.154 

8.299 

7.878 

7.705 

8.017 

8.249 

7.779 

6.717 

6.959 

7.354 

Results: 

average 

displacement  = 

9.187 

x  10-12  units 

calculated  value 

for  E 

=  2.177 

x  1011  Pascals 

=  31.74 

x  106  Psia 
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Test 


name:  PORE5-9 


Unit  length  =  1  micrometer  (10-*  m) 

•  80  cells,  4x4x5,  each  with  dimension  20x20x20  units 

Eo  =  46.7  x  10*  psia,  3.22  x  10u  Pascals 
22.03  %  porosity 

Pore  size:  cubes,  10  units  on  each  side 
Simulated  material:  Porous  Aluminium  Nitride  (AIN) 
Pores  follow  a  poison  distribution: 

®  Pores  per  cell:  01234567 

Number  of  cells:  16  19  23  14  6  2 


8 


Data: 

listed  below  are  the  displacements  of  the  25  node  points  on  the  top  of 
the  model  (each  x  10  *u  units). 

•  i  X  ->  Y 


1.102 

1.072 

1.118 

1.429 

0.996 

1.022 

1.083 

1.253 

1.432 

0.895 

1.033 

1.133 

1.137 

1.208 

1.083 

1.133 

1.084 

1.098 

1.170 

1.148 

0.988 

1.050 

1.128 

•  Results: 

average  displacement  =  10.844  x  10*12  units 
calculated  value  for  E  =  1.844  x  10u  Pascals 

=  26.748  x  10*  Psia 
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Test 


PORES- 1 

Unit  length  =  1  micrometer  (10-*  m) 

80  cells,  4x4x5,  each  with  dimension  20x20x20  units 
Eo  =  46.7  x  106  psia,  3.22  x  10u  Pascals 
3.574%  porosity 

Pore  size:  cubes,  5  units  on  each  side 
Simulated  material:  Porous  Aluminium  Nitride  (AIN) 

Pores  follow  a  poison  distribution: 

Pores  per  cell:  012345678 
Number  of  cells:  9  14  23  19  10  4  1 

Data: 

listed  below  are  the  displacements  of  the  25  node  points  on  the  top  of 
the  model  (each  x  10  ~n  units), 
i  X  ->  Y 


6.720 

6.731 

6.678 

6.573 

6.567 

6.718 

6.711 

6.657 

6.658 

6.731 

6.702 

6.730 

6.706 

6.741 

6.864 

6.729 

6.793 

6.767 

6.712 

6.743 

6.784 

6.807 

6.740 

6.638 

6.605 

Results: 

average  displacement  =  6.712  x  10-u  units 
calculated  value  for  E  =  2.980  x  10u  Pascals 

=  43.214  x  10«  Psia 


9  Test  name:  PORE6-2 

Unit  length  =  1  micrometer  (10-6  m) 

80  cells,  4x4x5,  each  with  dimension  20x20x20  units 
Eo  =  46.7  x  106  psia,  3.22  x  10u  Pascals 
5.41%  porosity 

Pore  size:  cubes,  5  units  on  each  side 
9  Simulated  material:  Porous  Aluminium  Nitride  (AIN) 

Pores  follow  a  poison  distribution: 

Pores  per  cell:  012345678 
Number  of  cells:  2  6  15  20  17  11  6  2  1 

Data: 

•  listed  below  are  the  displacements  of  the  25  node  points  on  the  top  of 

the  model  (each  x  10  _12  units), 
i  X  ->  Y 


7.162 

7.205 

7.148 

7.135 

7.241 

7.156 

7.128 

7.104 

7.153 

7.173 

7.035 

7.033 

7.095 

7.128 

7.094 

6.943 

6.984 

6.972 

6.947 

7.000 

6.965 

6.909 

6.765 

6.777 

6.671 

Results: 

average  displacement  =  7.045  x  10"12  units 
calculated  value  for  E  =  2.839  x  10u  Pascals 

=  41.17  x  106  Psia 
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Test 


name:  PORE6-3 

Unit  length  =  1  micrometer  (10-6  m) 

80  cells,  4x4x5,  each  with  dimension  20x20x20  units 
Eo  =  46.7  x  10®  psia,  3.22  x  10u  Pascals 
6.97  X  porosity 

Pore  size:  cubes,  5  units  on  each  side 
Simulated  material:  Porous  Aluminium  Nitride  (AIN) 

Pores  follow  a  poison  distribution: 

•  Pores  per  cell:  0123456789  10 

Number  of  cells:  1  2  8  15  17  15  11  6  3  1  1 


Data: 

listed  below  are  the  displacements  of  the  25  node  points  on  the  top  of 
the  model  (each  x  10  ~n  units), 
i  X  ->  Y 


7.163 

7.205 

7.157 

7.161 

7.275 

7.183 

7.221 

7.212 

7.228 

7.259 

7.213 

7.221 

7.253 

7.255 

7.235 

7.180 

7.204 

7.210 

7.193 

7.196 

7.142 

7.153 

7.138 

7.139 

7.132 

%  Results: 

average  displacement  =  7.197  x  10*12  units 
calculated  value  for  E  =  2.779  x  10u  Pascals 

=  40.303  x  10*  Psia 


Test  naae:  PORE6-4 

Unit  length  ~  1  micrometer  (10**  m) 

80  cells,  4x4x5,  each  with  dimension  20x20x20  units 
Eo  =  46.7  x  10*  psia,  3.22  x  10u  Pascals 
8.46  X  porosity 

Pore  size:  cubes,  5  units  on  each  side 
Simulated  material:  Porous  Aluminium  Nitride  (AIN) 

Pores  follow  a  poison  distribution: 

Pores  per  cell:  0123456789  10 

Number  of  cells:  0  1  4  9  14  15  14  10  7  4  2 

Data: 

listed  below  are  the  displacements  of  the  25  node  points  on  the  top  of 
the  model  (each  x  10  ~n  units). 

1  X  ->  Y 


7.322 

7.370 

7.376 

7.315 

7.289 

7.329 

7.348 

7.354 

7.334 

7.287 

7.349 

7.328 

7.328 

7.305 

7.273 

7.304 

7.304 

7.292 

7.257 

7.211 

7.255 

7.249 

7.223 

7.209 

7.184 

Results: 

average  displacement  =  7.296  x  10*12  units 
calculated  value  for  E  =  2.741  x  10u  Pascals 

=  39.757  x  10*  Psia 
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Teat 


name:  PORE&-5 

Unit  length  =  1  micrometer  (10-®  m) 

80  cells,  4x4x5,  each  with  dimension  20x20x20  units 
Eo  =  46.7  x  10®  psia,  3.22  x  1011  Pascals 
10.61  %  porosity 

Pore  size:  cubes,  5  units  on  each  side 
Simulated  material:  Porous  Aluminium  Nitride  (AIN) 

Pores  follow  a  poison  distribution: 

•  Pores  per  cell:  0123456789  10 

Number  of  cells:  0  0  1  3  7  11  13  14  13  11  7 

Data: 

listed  below  are  the  displacements  of  the  25  node  points  on  the  top  of 
the  model  (each  x  10  _l2  units). 

•  i  X  ->  Y 


7.492 

7.585 

7.611 

7.536 

7.497 

7.539 

7.560 

7.573 

7.547 

7.484 

7.524 

7.537 

7.538 

7.503 

7.472 

7.414 

7.451 

7.469 

7.458 

7.447 

7.332 

7.337 

7.344 

7.405 

7.450 

#  Results: 

average  displacement  =  7.484  x  10-12  units 
calculated  value  for  E  *  2.672  x  10u  Pascals 

=  38.757  x  10®  Psia 


Teat  name:  PORE6-6 

Unit  length  *  1  micrometer  (10*®  m) 

80  cells,  4x4x5,  each  with  dimension  20x20x20  units 
Eo  =  46.7  x  10®  psia,  3.22  x  10u  Pascals 
12.46  X  porosity 

Pore  size:  cubes,  5  units  on  each  side 
Simulated  material:  Porous  Aluminium  Nitride  (AIN) 

Pores  follow  a  poison  distribution: 

Pores  per  cell:  0123456789  10 

Number  of  cells:  0  0  0  1  2  5  8  12  16  18  18 


Data: 

listed  below  are  the  displacements  of  the  25  node  points  on  the  top  of 
the  model  (each  x  10  _12  units), 
i  X  ->  Y 


7.446 

7.498 

7.499 

7.446 

7.403 

7.538 

7.536 

7.517 

7.465 

7.403 

7.575 

7.561 

7.519 

7.469 

7.432 

7.535 

7.538 

7.512 

7.468 

7.444 

7.510 

7.493 

7.484 

7.456 

7.416 

Results: 

average 

displacement  = 

7.487 

x  10-12  units 

calculated  value  for  E 

=  2.671 

x  10u  Pascals 

=  38.744  x  10®  Psia 
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Test 


none:  PORE7-1 

Unit  length  =  1  micrometer  (10*®  m) 

80  cells,  4x4x5,  each  with  dimension  12x12x12  units 
Eo  =  46.7  x  106  psia,  3.22  x  10u  Pascals 
5.43  %  porosity 

Pore  size:  cubes,  5  units  on  each  side 
Simulated  material:  Porous  Aluminium  Nitride  (AIN) 

Pores  follow  a  poison  distribution: 

w  Pores  per  cell:  012345678 

Number  of  cells:  41  23  12  3  1 

Data: 

listed  below  are  the  displacements  of  the  25  node  points  on  the  top  of 
the  model  (each  x  10  *11  units). 

•  1  X  ->  Y 


1.254 

1.205 

1.184 

1.190 

1.183 

1.202 

1.171 

1.153 

1.159 

1.164 

1.124 

1.123 

1.119 

1.139 

1.169 

1.079 

1.092 

1.111 

1.133 

1.153 

1.075 

1.087 

1.114 

1.118 

1.107 

•  Results: 

average  displacement  =  11.44  x  10* 12  units 
calculated  value  for  E  =  2.914  x  10u  Pascals 

=  42.267  x  10®  Psia 


Test  name:  PORE7-2 

Unit  length  =  1  micrometer  ( 10*®  m) 

80  cells,  4x4x5,  each  with  dimension  12x12x12  units 
Eo  =  46.7  x  10®  psia,  3.22  x  10u  Pascals 
8.41  %  porosity 

Pore  size:  cubes,  5  units  on  each  side 
Simulated  material:  Porous  Aluminium  Nitride  (AIN) 

Pores  follow  a  poison  distribution: 

Pores  per  cell:  012345678 
Number  of  cells:  28  23  29  8  2 

Data: 

listed  below  are  the  displacements  of  the  25  node  points  on  the  top  of 
the  model  (each  x  10  *u  units), 
i  X  ->  Y 


1.334 

1.325 

1.292 

1.252 

1.229 

1.269 

1.256 

1.244 

1.246 

1.246 

1.173 

1.221 

1.238 

1.231 

1.269 

1.115 

1.199 

1.230 

1.228 

1.251 

1.120 

1.141 

1.211 

1.243 

1.222 

Results; 

average  displacement  =  12.31  x  10*12  units 
calculated  value  for  E  =  2.708  x  10u  Pascals 

=  39.277  x  10®  Psia 


Test  oame:  PORE7-3 

Unit  length  =  1  micrometer  (10~®  m) 

•  80  cells,  4x4x5,  each  with  dimension  12x12x12  units 

Eo  =  46.7  x  10®  psia,  3.22  x  1011  Pascals 
11.75  X  porosity 

Pore  size:  cubes,  5  units  on  each  side 
Simulated  material:  Porous  Aluminium  Nitride  (AIN) 
Pores  follow  a  poison  distribution: 

®  Pores  per  cell:  012345678 

Number  of  cells:  18  20  23  13  5  1 


Data: 

listed  below  are  the  displacements  of  the  25  node  points  on  the  top  of 
the  model  (each  x  10  -u  units). 

•  i  X  ->  Y 


1.509 

1.590 

1.587 

1.421 

1.350 

1.413 

1.435 

1.402 

1.362 

1.323 

1.288 

1.328 

1.339 

1.306 

1.307 

1.287 

1.365 

1.339 

1.331 

1.397 

• 

1.554 

Results: 

1.349 

1.337 

1.441 

1.468 

average 

displacement 

=  13.93 

x  10"12  units 

calculated  value  for  E  =  2.394  x  10u  Pascals 

=  34.72  x  10®  Psia 


•  Test  name:  PORE7-4 

Unit  length  =  1  micrometer  (10-®  m) 

80  cells,  4x4x5,  each  with  dimension  12x12x12  units 
Eo  =  46.7  x  10®  psia,  3.22  x  10u  Pascals 
14.65  %  porosity 

Pore  size:  cubes,  5  units  on  each  side 

•  Simulated  material:  Porous  Aluminium  Nitride  (AIN) 

Pores  follow  a  poison  distribution: 

Pores  per  cell:  012345678 
Number  of  cells:  12  17  24  16  7  3  1 

Data: 

•  listed  below  are  the  displacements  of  the  25  node  points  on  the  top  of 
the  model  (each  x  10  _u  units). 

IX  ->  Y 


1.578 

1.640 

1.628 

1.446 

1.359 

1.538 

1.533 

1.481 

1.435 

1.391 

1.493 

1.484 

1.448 

1.400 

1.391 

1.548 

1.561 

1.472 

1.425 

1.465 

1.828 

1.559 

1.487 

1.540 

1.530 

Results: 

average  displacement  =  15.06  x  10-12  units 
calculated  value  for  E  =  2.214  x  10u  Pascals 

=  32.11  x  10®  Psia 
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Test  Dane:  PORE7-5 

Unit  length  =  1  micrometer  (10-6  m) 

80  cells,  4x4x5,  each  with  dimension  12x12x12  units 
Eo  =  46.7  x  106  psia,  3.22  x  1011  Pascals 
18.17  X  porosity 

Pore  size:  cubes,  5  units  on  each  side 
Simulated  material:  Porous  Aluminium  Nitride  (AIN) 
Pores  follow  a  poison  distribution: 


Pores  per 

cell: 

0  1 

2  3 

4 

5  6  7 

8 

Number  of 

cells: 

7  12 

22  20 

12 

5  2 

Data: 

listed 

below 

are  the 

displacements  of 

the  25  node 

points  on  the  top  of 

the  model 

(each  x 

10  *u  units). 

1  X  -> 

Y 

1.967 

1.925 

1.730 

1.615 

1.771 

1.856 

1.777 

1.649 

1.598 

1.555 

1.872 

1.762 

1.687 

1.585 

1.464 

1.947 

1.830 

1.702 

1.608 

1.551 

1.999 

1.734 

1.646 

1.674 

1.644 

Results: 

average  displacement 

=  17.259 

x  10*12 

units 

calculated  value  for  E  s  1.932  x  10u  Pascals 

=  28.02  x  10*  Psia 


Test  name:  PORE7-6 

Unit  length  =  1  micrometer  (10-6  m) 

80  cells,  4x4x5,  each  with  dimension  12x12x12  units 
Eo  =  46.7  x  106  psia,  3.22  x  1011  Pascals 
20.44  X  porosity 

Pore  size:  cubes,  5  units  on  each  side 
Simulated  material:  Porous  Aluminium  Nitride  (AIN) 
Pores  follow  a  poison  distribution: 


Pores  per 

cell: 

0  1 

2  3 

4 

5  6 

7 

8 

Number  of 

cells: 

5  10 

20  20 

14 

7  3 

1 

Data: 

listed 

below  ; 

are  the  displacements  of 

the  25 

node 

points  on  the  top  of 

the  model 

(each  x 

10  *u  units). 

l  X  -> 

Y 

1.976 

1.933 

1.733 

1.619 

1.771 

1.904 

1.826 

1.695 

1.683 

1.788 

1.955 

1.848 

1.775 

1.728 

1.621 

2.057 

1.949 

1.831 

1.742 

1.691 

2.145 

1.888 

1.812 

1.852 

1.831 

Results: 

average  displacement 

=  18.263 

x  10*12 

units 

calculated  value  for  E  =  1.826  x  10u  Pascals 

=  26.48  x  106  Psia 
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Test 


P0RE8-1 

Unit  length  =  1  micrometer  (1(H  m) 

•  80  cells,  4x4x5,  each  with  dimension  20x20x20  units 
Eo  =  46.7  x  10*  psia,  3.22  x  1011  Pascals 

5.2734  X  porosity 

Pore  size:  cubes,  15  units  on  each  side 
Simulated  material:  Porous  Aluminium  Nitride  (AIN) 

Pores  follow  a  poison  distribution: 

Pores  per  cell:  012345678 
Number  of  cells:  71  8  1 

Data: 

listed  below  are  the  displacements  of  the  25  node  points  on  the  top  of 
the  model  (each  x  10  ~n  units). 

•  i  X  ->  Y 


6.528 

6.488 

6.428 

6.338 

6.197 

6.649 

6.677 

6.672 

6.512 

6.304 

7.232 

6.964 

6.732 

6.576 

6.376 

7.372 

6.985 

6.868 

6.754 

6.393 

6.943 

6.823 

7.184 

7.074 

6.419 

Results: 

average  displacement  =  6.700  x  10"12  units 
calculated  value  for  E  =  2.985  x  10u  Pascals 

=  43.295  x  10«  Psia 


•  Test  name:  PORE8-2 

Unit  length  r  1  micrometer  (10~®  m) 

80  cells,  4x4x5,  each  with  dimension  20x20x20  units 
Eo  =  46,7  x  106  psia,  3.22  x  10u  Pascals 
7.910  X  porosity 

Pore  size:  cubes,  15  units  on  each  side 

•  Simulated  material:  Porous  Aluminium  Nitride  (AIN) 

Pores  follow  a  poison  distribution: 

Pores  per  cell:  012345678 
Number  of  cells:  67  11  2 

Data: 

•  listed  below  are  the  displacements  of  the  25  node  points  on  the  top  of 
the  model  (each  x  10  *12  units). 

i  X  ->  Y 


6.412 

6.427 

6.409 

6.348 

6.253 

6.654 

6.749 

6.810 

6.717 

6.580 

7.345 

7.153 

7.015 

6.995 

6.928 

7.571 

7.275 

7.264 

7.291 

7.079 

7.230 

7.210 

7.667 

7.619 

7.073 

Results: 

average  displacement  =  6.963  x  10-12  units 
calculated  value  for  E  =  2.872  x  1011  Pascals 

=  41.657  x  106  Psia 
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Test  nase:  PORB8-3 

Unit  length  =  1  micrometer  (10-6  m) 

80  cells,  4x4x5,  each  with  dimension  20x20x20  units 
Eo  =  46.7  x  106  psia,  3.22  x  1011  Pascals 
11.074  %  porosity 

Pore  size:  cubes,  15  units  on  each  side 
Simulated  material:  Porous  Aluminium  Nitride  (AIN) 
Pores  follow  a  poison  distribution: 

Pores  per  cell:  012345678 
Number  of  cells:  62  15  3 


Data: 

listed  below  are  the  displacements  of  the  25  node  points  on  the  top  of 
the  model  (each  x  10  *12  units). 

IX  ->  Y 


7.587 

7.376 

7.158 

7.093 

7.069 

7.355 

7.340 

7.318 

7.194 

7.039 

7.589 

7.392 

7.198 

7.518 

8.304 

7.543 

7.243 

7.171 

7.562 

8.052 

6.968 

6.921 

7.375 

7.162 

6.470 

Results: 

average  displacement  =  7.320  x  10' 12  units 
calculated  value  for  E  =  2.732  x  10u  Pascals 

=  39.626  x  10«  Psia 


•  Test  name:  PORE8-4 

Unit  length  =  1  micrometer  (10*6  m) 

80  cells,  4x4x5,  each  with  dimension  20x20x20  units 
Eo  =  46.7  x  106  psia,  3.22  x  10u  Pascals 
14.238  X  porosity 

Pore  size:  cubes,  15  units  on  each  side 

•  Simulated  material:  Porous  Aluminium  Nitride  (AIN) 

Pores  follow  a  poison  distribution: 

Pores  per  cell:  012345678 
Number  of  cells:  59  17  5 

Data: 

®  listed  below  are  the  displacements  of  the  25  node  points  on  the  top  of 

the  model  (each  x  10  units). 

IX  ->  Y 


8.076 

7.886 

8.808 

8.426 

6,917 

7.728 

7.751 

8.017 

7.745 

7.112 

7.886 

7.687 

7.427 

7.756 

8.607 

7.783 

7.452 

7.388 

7.857 

8.438 

7.146 

7.093 

7.552 

7.366 

6.702 

Results: 

average  displacement  =  7.704  x  10'12  units 
calculated  value  for  E  =  2.596  x  10u  Pascals 

=  37.649  x  10«  Psia 
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Test 


name:  PORE8-5 

Unit  length  =  1  micrometer  (10-6  m) 

80  cells,  4x4x5,  each  with  dimension  20x20x20  units 
Eo  =  46.7  x  106  psia,  3.22  x  1011  Pascals 
18.984  %  porosity 

Pore  size:  cubes,  15  units  on  each  side 
Simulated  material:  Porous  Aluminium  Nitride  (AIN) 
Pores  follow  a  poison  distribution: 

Pores  per  cell:  012345678 
Number  of  cells:  53  20  8 


Data: 

listed  below  are 
the  model  (each  x  10 
IX  ->  Y 


the  displacements 
~12  units). 


of  the  25  node  points  on  the  top  of 


7.920 

7.869 

9.062 

10.081 

9.847 

7.606 

8.100 

8.711 

8.802 

8.936 

7.817 

8.293 

8.212 

8.453 

9.513 

7.790 

7.837 

8.008 

8.621 

9.403 

7.104 

7.362 

8.067 

8.101 

7.648 

Its: 

average 

displacement 

=  8.3667 

x  10" 12  units 

calculated  value  for  E 


=  2.3904  x  10u 
=  34.669  x  10« 


Pascals 

Psia 
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A  FINITE  ELEMENT  ANALYSIS  OF  POROSITY  EFFECTS  ON  MATERIALS 
This  paper  examines  the  effects  of  porosity  using  a  three  dimensional  finite 
element  njodel.  Three  previously  used  analytical  equations  (linear,  empirical 
exponential,  and  Hassin’s  semi-empirical  equation)  are  presented.  These 
analytical  models  all  require  that  the  effects  of  porosity  on  materv_l  properties 
be  determined  experimentally  for  several  cases,  then  extrapolated  i.r  other 
cases.  This  limits  accurate  prediction  before  a  material  is  produced.  In  this 
paper  a  finite  element  model  is  developed  that  can  predict  the  material 
modulus  of  elasticity,  Young’s  modulus,  using  the  limited  information  from  one 
material  sample.  The  model  is  shown  to  correspond  with  several  known  test 
results.  The  same  finite  element  model  is  used  to  demonstrate  the  effects  of 
pore  size  and  distribution  on  material  behavior. 


